AO>AOM  217 


UNCLASSIFIED 


aerospace  medical  research  lab  WRI6HT-PATTERS0N  AFB  OHIO  F/0  6/2 

aerodynamic  forces  exerted  on  an  articulated  boot  subjected  to  -^TC(U) 

dec  76  D J SCHNECK 

AMRL-TR-76-109  NL 


1 W'l 

AO 

A0442I7 

■ 

1 

— 

j 

J 

i 

L 

- 

, '■ 

1 

j 

[ 

L 

: 1 

1 

1 

[ 

1 

i 

f 

END 

DATE 

niNED 

10-77 

MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURIAU  Of  STANDARDS  1963*.^ 


g AERODYNAMIC  FORCES  EXERTED  ON  AN  ARTICULATED 
BODY  SUBJECTED  TO  WINDBLAST 


AEROSPACE  MEDICAL  RESEARCH  LABORATORY 


CJaerospace  medical  research  laboratory 
Cia  ^Saerospace  medical  division 
£Z3air  force  systems  command 


WRIGHT-PATTERSON  AIR  FORCE  BASE.  OHIO  45433 


NOTICES 


When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than  a definitely  relau 
Government  procurement  operation,  the  Government  thereby  incurs  no  responsibility  nor  any  obligation  ah„ 
soever,  and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the  said  drawings, 
specifications,  or  other  data,  is  not  to  be  regarded  by  implication  or  otherwise,  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to  manufacture,  use.  or  sell  any  patented 
invention  that  may  in  any  way  be  related  thereto. 


Please  do  not  request  copies  of  this  report  from  Aerospace  Medical  Research  Laboratory.  Additional  copies  may  be 
purchased  Irom; 

National  Technical  Information  Service 
5285  Port  Royal  Road 
Springfield,  Virginia  22161 

Federal  Government  agencies  and  their  contractors  registered  with  Defense  Documentation  Center  should  direct 
requests  for  copies  of  this  report  to: 

Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 


TECHNICAL  REVIEW  AND  APPROVAL 

AMRL-TR-76-109 


This  report  has  been  reviewed  by  the  Information  Office  (Ol)  and  is  releasable  to  the  National  Technical  Information 
Service  (NTIS).  At  NTIS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


Jf. — 

HENNIN^^VON  GIERKE 
Director 


Biodynamics  and  Bionics  Division 
Aerospace  Medical  Research  Laboratory 


1 


AID 


rO«CC  - 4 FEHWIAPY  77  . 100 


i 

( 

! 


security  classification  of  This  page  (Wh»n  Dmfa  Entered} 


REPORT  DOCUMENTATION  PAGE 


' M-)  AMRL-TR-76-lg(9j 


5ArP  READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

2 GOVT  ACCESSION  NO.  ,-^i^eii.lPieNT'S  CATALO^  NUMBER 

i'O  - 


14.  Title  (and  Subtltta) 


>ERODYNAMIC  JORGES  EXERTED  ON  AN, ARTICULATED 

W-  ^ 

_BODY  SUBJECTED  TO  WIND  BLAST  , 

■ — ■ I 

7'  AuTtlg.BIii  . 

Daniel  J./Schneck/* 


^ irrTCPyRY'a'pgffioD  eo've~RfB" 
Final  ReportI  June  14- 
August  20,  1976 

6.  PtRFOJiMlMa  ORU.  REPORT  NUMBER 

.yy 


I ^enTRACT  OR  grant  NUMBEJy^; 


19.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Aerospace  Medical  Research  Laboratory,  Aerospace 
Medical  Division,  Air  Force  Systems  Command, 
Wriaht-Patterson  Air  Force  Base,  Ohio  45433 


111.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


y yi.  / 


10.  program  ELEMENT./PROJECT.  task 
AREA  a UNIT  ^UMBERS 

61102fJ2312-1v3-11 


14  MONITORING  AGENCY  NAME  a AODRES^f/f  di/ferent  from  ControWnH  Of/ice)  IS-  SSCUFlT  Y t L ASS.  repotr^ 

^ Unclassified 

~TSa  DECLASSIFICATION  OOWNGRADy 
SCHEDULE  f 


Ma.  DISTRIBUTION  ST ATEMENT  (of  this  Peprirt) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (o(  the  abstract  er\tered  In  Block  20.  U different  from  Report) 


<'■  yJ 

V- 


18  supplementary  notes 

* The  author  is  an  Associate  Professor,  Virginia  Polytechnic  Institute  and  State 
University,  Blacksburg,  Virginia  24061.  He  was  in  residence  in  AMRL  as  a 
Summer  Faculty  Fellow  under  a research  associateship  from  the  USAF/ASEE. 


I 19  key  WORDS  I'Conflnu©  on  reverse  side  if  necessary  and  identify  by  block  number) 


Biological  and  Medical  Sciences 
Biodynamic  Modeling 
Aerodynamic  Forces 
Inviscid  Flow  Theory 
jHuman  Body  Model 


Windblast  Injuries 
Flail  Injury 
High-Speed  Ejection 
Flow  Separation 
Potential  Flow  Theory 


^ ABSTRACT  fConffnue  on  reverse  aide  If  necessary  and  Identify  by  block  number) 

A potential  flow  solution  is  presented  for  estimating  the  pressure  distribution 
around  the  forearm  of  a human  body  subjected  to  windblast.  The  forearm  is  ex- 
amined in  three  positions:  resting  and  pressing  against  an  arm  rest;  resting, 
but  not  pressing  against  an  arm  rest;  and  not  resting  at  all  against  any  sur- 
face. Results  show  that  a high-speed  wind  stream  approaching  the  limb  at  some 
finite  angle  of  attack  has  a tendency  to  dislodge  the  forearm  from  a surface 
with  which  it  is  in  contact.  This  is  due  to  the  generation  of  stagnation  points 
in  the  flow  which  lead  to  adverse  pressure  gradients  as  high  as  six  times  the  ^ 

DD  . 1473  EDITION  OF  I NOV  65  IS  OBSOLETE 


security  CLASSIFICATION  OF  THIS  PAGE  C»»>ien  Entered 


r ^/'  r 


security  classification  of  this  PAQEfmi»n  D*l«  Enl«r«d; 


/20/  free-stream  dynamic  pressure.  Moreover,  when  the  inviscid  analysis  is 
corrected  for  the  effects  of  flow  separation,  it  is  possible  to  predict  the 
presence  of  a pressure  drag  which  acts  to  throw  the  forearm  outward,  away  from 
the  thorax.  Botli  of  these  effects  increase  with  angle  of  attack  and  they  are 
mildly  dependent  on  the  taper  of  the  forearm  geometry.  It  is  anticipated  that  a 
straightforward  extension  of  the  theory  to  other  limbs  of  the  body  will  shed 
some  light  on  the  general  problem  of  flail  injury  occurring  during  high-speed 
ej  ect ions . 
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SUMMARY 
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Problem 

The  specific  problem  was  to  formulate  a mathematical  approach  which 
will  assist  in  predicting  the  aerodynamic  loading  to  which  a pilot  is 
exposed  during  ejection  from  a high-perf omance  aircraft.  Such  aerody- 
namic data  is  intended  to  serve  as  input  information  to  the  Aerospace 
Medical  Research  Laboratory's  Articulated  Total  Body  Model  (ATBM) , which 
is  a further  development  and  refinement  of  the  Calspan  model.  This  model, 
subjected  to  suitably  defined  aerodynamic  loading,  can  be  used  effectively 
to  study  injuries  which  take  place  as  a result  of  limb-flailing  - that  is, 
the  tendency  of  the  arms  and  legs  to  swing  freely  and  out  of  control  due  to 
windblast  forces  that  exceed  musculoskeletal  resistance.  The  results  of  the 
following  study  will  also  guide  in  planning  significant  wind-tunnel  experi- 
ments which  will  provide  important  data  for  the  design  of  high-speed  ejec- 
tion mechanisms.  The  mathematical  model  as  developed  in  this  report  is 
intended  to  be  ultim.ateiy  correlated  with  wind-tunnel  experiments  in  order 
to  provide  a theoretical  framework  from  which  some  physical  insight  may  be 
gained  into  the  general  problem  of  limb-flailing. 

Approach 

The  body  me  lei  used  tor  the  analysis  of  windblast  forces  in  this  report 
consists  of  a IS-linkasC  system  composed  of  spherical,  circular-cylindrical, 
truncated-conical  and  flat-plate  segments  attached  at  a series  of  arbitrary 
joints.  Each  of  the  segments  can  be  described  geometrically  from  available 
anthropometric  data.  Its  orientation  in  space  can,  in  turn,  be  defined  by 
specifying  the  direction  cosines  of  a coordinate  system  fixed  to  the  limb, 
relative  to  an  inertial  coordinate  system  fixed  in  space.  To  illustrate  the 
feasibility  of  the  approach,  attention  has  been  focused  specifically  on  the 
forces  acting  on  the  forearm  portion  of  the  right  arm.  A free  stream  is 
assumed  to  approach  the  limb  at  some  angle  of  attack  relative  to  the  fore- 
arm centerline.  The  latter  is  examined  in  three  specific  positions:  rest- 

ing and  pressing  against  an  armrest;  resting,  but  not  pressing  against  an 
armrest;  and  not  resting  at  all  against  any  surface.  For  practical  purposes, 
the  stream  was  considered  to  be  incompressible  and  inviscid,  so  that  a stream 
function  and  velocity  potential  could  be  defined  for  each  of  the  three  situ- 
ations described.  In  addition,  turbulence  and  three-dimensional  effects 
have  been  temporarily  neglected. 

Results 

Solution  of  the  conservation  of  mass  and  momentum  equations  allovjed  the 
pressure  distribution  around  the  forearm  to  be  calculated.  For  the  case  of 
the  forearm  not  resting  against  any  surface,  it  was  demonstrated  that  pres- 
sure maxima  occur  symmetrically  at  the  front  and  back  of  the  limb,  corres- 
ponding to  points  of  flow  stagnation.  These  maxima  decreased  with  angle  of 
attack.  It  was  also  found  that  pressure  minima  occur  symmetrically  on  the 
top  and  bottom  of  the  forearm  when  it  is  at  right  angles  to  the  incident 
stream,  but  the  points  themselves  move  leeward  with  decreasing  angle  of 
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attack  at  a rate  mildly  dependent  on  the  hall-angle  ol  tlie  t rnnca t ed-enne 
torearm  font igurat ion. 

For  the  case  ot  the  forearm  resting  against  a solid  surface,  it  was 
demonstrated  tiiat  pressure  maxima  occur  at  the  point  of  contact  In'tween  the 

limb  and  the  flat  restraining  plate,  and  that  there  is  a high  pressure  region  1 

spanning  nearly  160  degrees  of  arc  along  the  lower  half  of  the  circular  cross  j 

section.  Moreover,  there  is  a relatively  deep  low  pressure  region  spanning  j 

nearly  200  degrees  of  arc  along  the  top  of  the  forearm,  such  that  there  I 

results  a rather  substantial  net  force  that  acts  to  lift  the  cone  (or  dis-  1 

lodge  the  limb)  from  the  surface  with  which  it  is  in  contact.  It  was  also 

found  that  pressure  minima  occur  near  the  top  of  the  limb  and  move  leeward 
with  decreasing  angle  of  attack  at  a rate  mildly  dependent  on  the  cone  half- 
angle. The  net  force  tending  to  cause  limb  dislodgement  was  computed  to  be  ’ 

some  six  times  the  free  stream  dynamic  pressure. 

I 

When  the  inviscid  analysis  was  corrected  for  the  effects  of  flow  separa-  1 

tion,  it  was  possible  to  argue  qualitatively  for  the  existence  of  a profile, 

form  or  pressure  drag  which  acts  to  throw  the  forearm  outward,  away  from  the  1 

thorax.  The  magnitude  of  this  tendency  is  yet  to  be  quantified,  but  it  is  in  | 

agreement  with  wind  tunnel  results  which  have  shown  a distinct  upward  and  ij 

outward  motion  of  the  forearm  when  subjected  to  conditions  similar  to  those  J 

simulated  here.  j 


Conclusions  < 

{ 

It  has  been  shown  that  one  can  approach  the  general  problem  of  limb-  ! 

flailing  during  ejection  by  constructing  a simplified  model  from  which  the  . 

aerodynamic  forces  exerted  on  an  articulated  body  subjected  to  windblast  may  ; 

be  reasonably  estimated.  This  model  has  clearly  demonstrated  that  flow  j 

separation,  together  with  the  appearance  of  stagnation  points  in  the  fluid  , 

present  themselves  as  first  order  effects  responsible  for  the  generation  of  1 

limb-dislodging  forces  that  result  in  flail  injuries.  The  former  gives  rise  t 

to  a pressure  gradient  that  tends  to  drag  the  limb  along  with  the  air  flow,  ; 

and  the  latter  results  in  a pressure  gradient  that  tends  to  separate  the  limb  ■ 

from  a surface  with  which  it  is  in  contact.  Both  of  these  effects  are  func-  t 

tions  of  the  square  of  the  velocity  of  the  incident  stream,  and  so  they 

become  rather  significant  at  speeds  near  sonic.  The  model  developed  here  can  'i 

be  progressively  refined,  made  more  sophisticated,  and  stream.llned  to  provide  , 

important  data  for  both  experimental  and  numerical  studies  of  human  body  J 

dynamics.  ( 

i 

Recommendat ions 

In  the  early  form  developed  thus  far,  the  body  joints  themselves  have 
not  been  defined,  and  it  would  be  desirable  to  look  into  this  in  greater  de-  1 

tail  in  order  to  assess  the  significance  of  finite  length  segments,  end  ef- 
fects, three-dimensional  contours  and  interactions  wliich  take  place  between  I 

the  limbs  as  a result  of  tlieir  joint  connections.  Along  the  same  lines,  the 
velocity  approaching  any  limb  must  be  viewed  as  the  resultant  I'f  an  outer 
flow  that  has  undergone  specific  interactions  witli  other  body  segment^  before 
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it  actually  inipiiigcs  ui)on  the  limb,  and  this  is  an  area  that  has  not  yot 
been  explored. 


In  retining  the  tlieory,  one  might  also  back-track  to  examine  more 
closely  tile  consequences  of  other  assumptions  made  in  the  analysis.  For 
e.xample,  the  rougliness  of  the  surface  over  which  the  fluid  flows  — which 
may,  in  part,  be  due  to  the  type  of  clothing  being  worn  by  the  pilot  — may 
lead  to  certain  types  of  turbulent  beliavior  that  is  not  negligible. 
Similarly,  at  .Mach  numbers  exceeding  0.5,  compressibiJ ity  effects  should 
certainly  be  included  and  flow  separation  effects  must  be  accurately  quanti- 
f ied . 
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SKCTION  1 


INTRODUCTION 


It  lias  been  estimated  (Pavne,  1975)  that  nearly  half  of  all  ejections 
taking  place  under  combat  conditions  result  in  flail  injury  or  death  to  the 
pilot.  This  is  therefore  a problem  of  some  magnitude  and  it  must  be  dealt 
with  accordingly.  In  an  effort  to  do  so,  the  Aerospace  Medical  Research 
Laboratory  of  the  Air  Force  Systems  Command  is  pursuing  experimental  and 
analytic  investigations  into  the  dangers  of  ejecting  a human  body  into  a 
high-speed  wind  stream.  Some  preliminary  experimental  work  in  this  area  has 
already  been  reported  (Hawker  and  Euler,  1975,  Brinkley  and  Payne,  1973)  and 
more  wind  tunnel  tests  are  being  planned. 

The  analytic  studies  underway  at  AMRL  are  utilizing  a computerized 
Articulated  Total  Body  Model  (ATBM)  to  describe  the  dynamics  of  the  human 
body  under  a variety  of  situations.  The  ATBM  has  grown  out  of  the  original 
Calspan  Model  (Fleck,  Butler  and  Vogel,  1974,  Bartz,  1971,  Bartz  and  Butler, 
1972)  altered  to  include  improved  joint  formulation,  improved  harness  formu- 
lation and  the  addition  of  aerodynamic  forces  (Fleck  and  Butler,  1975).  The 
latter,  however,  have  not  been  clearly  defined  because  little  is  known  about 
the  aerodynamic  loading  to  which  a pilot  is  exposed  during  ejection  from  a 
high  performance  aircraft  at  near  sonic  speeds.  The  work  described  herein 
was  undertaken  to  provide  some  initial  mathematical  basis  for  analyzing  and 
estimating  such  loading.  The  results  are  intended  primarily,  and  ultimately 
to  serve  as  input  data  to  the  ATBM  Program.  However,  they  will  also  guide 
and  be  correlated  with  corresponding  wind  tunnel  experiments  to  provide  a 
theoretical  framework  from  which  some  physical  Insight  may  be  gained  re- 
garding certain  mechanisms  which  cause  the  limbs  to  swing  freely  and  out  of 
control,  i.e.,  to  "flail"  during  ejection  due  to  windblast  forces  that 
exceed  musculo-skeletal  resistance. 


SECTION  II 

THE  BODY  MODEL  USED  FOR  THE  ANALYSIS  OF  WINDBLAST  FORCES 


For  the  purposes  of  this  analysis,  the  human  organism  has  been  modeled 
as  a 15-linkage  system  composed  of  spherical,  circular-cylindrical,  trun- 
cated-conical and  flat-plate  segments  attached  through  a series  of  arbi- 
trary joints  (as  yet,  not  specifically  defined).*  A schematic  diagram  of 
this  model  is  shown  in  Figure  1,  where: 


* The  body  model  being  developed  here  for  the  study  of  windblast  forces  dif- 
fers from  the  ATBM  in  the  number  of  segments  employed  and  in  their  shape. 
The  ATBM  uses  ellipsoids  throughout  the  program  for  modeling  the  contact 
surfaces  of  the  body  and  other  curved  surfaces  such  as  interior  surfaces, 
etc . 


1 


Segment 


1 Head  (Sphere)  ; 

2 = Neck  (Circular  Cylinder); 

3,  A = Dp[)er  Arms  (Circular  Cylinders); 

5,6  = Forearms  (Truncated  Circular  Cones); 

7 = Hand  Gripping  Seat  (Fist,  Sphere); 

8 = Hand,  Fingers  Extended  (Flat  Plate); 

9 = Torso  (Circular  Cylinder,  or  Elliptic  Cylinder); 

10,11  = Thighs  (Truncated  Circular  Cones); 

12,13  = Shanks  (Truncated  Circular  Cones); 

14,15  = Feet  (Flat  Plates). 

Geometric  description  of  each  of  the  15  segments  is  made  possible  by  ex- 
isting anthropometric  data  (Garret  and  Kennedy,  1971,  Clauser,  et  al,  1972). 
Its  respective  orientation  in  space,  in  turn,  is  defined  by  specifying  the 
direction  cosines  of  the  axes  of  a coordinate  system  fixed  to  the  limb, 
relative  to  an  inertial  coordinate  system  fixed  in  space.  Consider,  for 
example,  the  x-y-z  coordinate  system  fixed  to  the  right  forearm  (segment  5) 
as  illustrated  in  Figure  2.  The  x-axis  is  chosen  to  be  the  center  line  of 
the  truncated  circular  cone;  the  y-axis  is  oriented  so  that  the  free  stream 
approaches  the  limb  parallel  to  the  x-y  plane,  and  the  z-axis  lies  perpen- 
dicular to  the  x-y  plane.  The  position  of  the  forearm  may  be  defined  at  ai>y 
time,  t,  by  locating  a unit  vector  along  the  x-axis  relative  to  the  inertial 


coordinate  system,  x'-y'-z' 


Thus,  if  Xj^'(t),  y^'(t)  and  z^'(t)  represent 


the  location  of  the  elbow  Joint,  and  X2'(t),  y2'(t)  and  z.^' (t)  represent  the 

location  of  the  wrist  joint  at  any  time,  t,  in  the  x'-y'-z'  coordinate  sys- 
tem, then  the  vector  from  A to  B along  the  x-axis  is  defined  as  follows: 


r^(t)  + AB  = r2(t) 


from  which, 

AB  = r2(t)  - ?^(t). 


AB  = [x^'(t)  - x^'(t)Ji  + [y2'(t)  - y^'CtlJj  + - z^'(r)Ik  (2) 

wliere  i,  j,  and  k represent  unit  vectors  along  the  x',  y'  and  z'  directions, 
respectively.  Now,  the  magnitude  of  the  vector  AB  is  the  square  root  of  the 
sum  of  tlie  sqviares  of  its  components,  i.e.. 


|AB|  = I 1x2'  (t)  - x^  ' (t) 


+ [y2'(>-)  - y^'Ct)]' 


+ lz/(t)  - z^'(t)]^} 


Tlius,  a unit  vector  along  tlie  x-axis  (the  forearm  centerline)  t.l  .tl.c  tc 
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tile  inertial  coordinate  system  is  given  simply  by  equation  (2)  divided  by 
equation  (3),  that  is,  the  vector  AB  divided  by  its  magnitude: 


\ = unit  vector  from  elbow  to  wrist 

X 


AB(t) 

|A^| 


(t). 


(4) 


The  components  of  this  unit  vector  are  the  direction  cosines  of  the  angle 
that  the  x-axis  makes,  respectively,  with  the  x',  y'  and  z'  axes.  In  a 
similar  manner,  one  may  go  on  to  completely  describe  the  spatial  position 

of  any  body  segment,  i,  at  any  time,  t,  by  defining  the  unit  vectors 

(i  = 1 to  15,  j = x, y, z) attached  to  each  limb  at  a specified  joint: 


Limb  position  (general)  = A^^^t). 


(5) 


Now,  if  a free  stream  approaches  the  forearm  with  some  constant  veloc- 


ity, 


U = L'.i  + U_j  + U_ k = U A„ 
1 2.  j O u 


(6) 


such  that: 


A = — = unit  vector  in  the  stream  direction, 
u L 


(7) 


then  the  "angle  of  attack"  between  the  oncoming  stream  and  the  forearm 
centerline  may  be  defined  to  be  a,  where: 


cos  a (t)  = ^jj(t)  • Ajj 


(8) 


Note,  from  Figure  2, 

stream  U can  be  made 
such  that: 


that  by  properly  orienting  the  y-axis,  the  incident 
to  lie  in  the  x-y  plane,  having  components  -U  and  -V 


tan  a = . (9) 

U is  then  the  velocity  component  directed  parallel  to  the  limb  centerline, 
and  V is  generally  referred  to  as  the  "cross-flow"  ve  ocity  component  in  a 
plane  perpendicular  to  the  x-axis. 

Continuing  to  concentrate,  in  this  initial  analysis,  on  the  right  fore- 
arm, we  introduce  a cylindrical  coordinate  system,  r,  6,  x as  shown  in 
Figure  2 and  define  the  cross-sectional  radius  of  the  forearm  at  any 
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location  x measured  from  the  elbow  as  R.  Then  it  is  clear  that  R = R(x) 
o 

and,  if  the  forearm  is  modeled  as  a truncated  right  circular  cone, 


tan  3 = = Cone  Half-Angle.  (10) 

dx 

Note  that,  since  R(x)  decreases  with  increasing  x the  way  the  coordinate 
system  has  been  set  up  in  the  limb,  tan  3 as  defined  by  equation  (10)  is 
negative.  Similarly,  tan  a as  defined  by  equation  (9)  is  positive,  since  V 
and  U are  both  assumed  to  be  negative.  Finally,  the  angle  6 is  considered 
positive  when  measured  counterclockwise  from  the  positive  y-axis.  As  men- 
tioned previously,  the  parameters  3 (cone  half-angle),  I (total  forearm 
length  from  elbow  to  wrist  = |AB|),  R(x)  and  others  pertinent  to  the  geo- 
metric description  of  the  limb  may  be  deduced  from  available  anthropometric 

-► 

data  for  adult  men  and  women.  The  free  stream  velocity,  U refers  specifi- 
cally to  the  velocity  approaching  the  forearm  "as  the  forearm  sees  it," 
i.e.,  including  the  effects  of  the  presence  of  a torso  and  other  wind-de- 
flecting disturbances.  It  is  thus  intimately  coupled  with  the  flow  regimes 
surrounding  the  other  body  segments  and  its  determination  must  necessarily 
involve  solving  an  entire  complex  system  of  equations  simultaneously.  This 
matter  shall  be  dealt  with  in  greater  detail  during  a subsequent  phase  of 
this  investigation.  For  the  present,  we  assume  only  the  simplest  case  - 

-V 

that  is,  that  U is  the  resultant  of  the  average  velocity  of  the  pilot  with 
respect  to  the  plane  during  the  ejection  phase  plus  the  velocity  of  the 
plane  with  respect  to  the  inertial  reference  frame  (e.g.,  the  ground): 


U = U -T  ./  1 + 

pilot/ plane 


U 


plane/ground 


(11) 


Some  additional  comments  relative  to  equation  (11)  shall  be  made  in  the 
Concluding  Remarks  (Section  VIII) . 


SECTION  III 

PRELIMINARY  ASSUMPTIONS 


In  the  most  general  sense,  this  problem  involves  the  turbulent  flow  of 
a viscous,  compressible  fluid  around  a complicated  three-dimensional  15- 
linkage  system  of  finite  bodies  of  revolution  attached  at  undefined  joints. 
Obviously,  one  must  begin  modestly  by  making  significant  simplifying  assump- 
tions if  any  degree  of  mathematical  tractability  is  to  be  sought  in  getting 
an  order  of  magnitude  estimate  of  the  aerodynamic  forces  with  which  we  are 
deal Ing. 

The  first  assumption  which  clearly  simplifies  matters  is  the  neglect 
of  turbulence,  i.e.,  the  flow  is  considered  to  be  laminar  except,  perhaps, 
in  a finite  region  on  the  downstream  side  of  the  bluff  body  segments  where 
flow  Separation  occurs.  The  laminar-flow  hypothesis  is  introduced  for  the 
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sake  of  convenience.  Its  validity  can  t'l-  i-valuated  by  comparing  analytic 
results  with  experimental  data,  which  will  he  done  as  such  data  becomes 
available . 

The  second  assumption  is  certainly  more  controversial,  but  has  led  to 
some  surprisingly  meaningful  results  when  applied  to  analogous  flow  patterns 
over  slender  inclined  bodies  of  revolution  (see,  for  example,  Allen  and 
Perkins,  1951,  Kelly,  1954  and  Mello,  1959).  That  is,  it  shall  be  assumed 
for  the  moment  that  the  effects  of  viscosity  in  the  fluid  may  be  linearly 
super-imposed  on  the  invlscid  flow  solutions,  such  that  the  latter  may  be 
determined  as  if  they  were  completely  independent  of  any  surface  viscous 
effects.  Looking  at  the  problem  in  greater  detail,  one  observes  that  this 
assumption  is  probably  not  as  bad  as  it  may  sound.  In  general,  during  the 
early  stages  of  flow  development  around  a blunt  body  of  revolution,  the 
vorticity  in  the  fluid  is  practically  confined  to  the  boundary  layer  and  to 
the  separated  wake  region  behind  the  body  (Goldstein,  1938,  page  61).  Now 
in  a very  high  Reynolds  number  flow  the  laminar  boundary  layer  can  be  ex- 
pected to  be  quite  thin,  and,  during  the  first  few  milliseconds  of  flow 
development,  thinner  still  because  not  much  time  has  elapsed  to  allow  the 
boundary  layer  to  grow.  It  thus  seems  likely  that  during  this  time  the 
wake  region  dominates  the  flow  patterns  to  a much  greater  extent  than  do  any 
effects  taking  place  in  the  growing  surface  boundary  layer  itself.  In  fact, 
even  in  fully-developed  flow  situations  over  bluff  bodies,  the  skin-friction 
drag  is  small  compared  with  the  form  drag,  or  pressure  drag,  due  to  flow 
separation  (Goldstein,  1938,  page  65).  Considering  that  the  problem  being 
dealt  with  here  is  one  of  very  high  Reynolds  number  (near  sonic  velocities 
of  air  the  Reynolds  number  is  of  order  10^) , and  that  most  flail  injuries 
take  place  within  the  first  3 milliseconds  of  the  ejection  process,  it 
would  seem  that,  for  practical  purposes,  the  flow  up  to  the  point  of  bound- 
ary-layer  separation  may  be  considered  to  be  inviscid  without  destroying  the 
essential  physics  of  the  problem. 

Going  one  step  further,  the  pressure  drag  due  to  flow  separation  may  be 
accounted  for  by  a super-position  of  invlscid  solutions  - if  the  boundary 
layer  on  the  surface  of  the  body  is  neglected  following  the  reasoning  of  the 
previous  paragraph.  Consider,  for  example,  the  case  of  a circular  cylinder 
placed  so  that  its  axis  is  at  a right  angle  to  an  oncoming  stream. 

Goldstein  (1938,  page  60)  points  out  that  when  the  flow  is  started  from 
rest,  separation  of  the  boundary  layer  begins  at  the  rear  of  the  cylinder 
after  a period  of  time  that  is  roughly  proportional  to  the  radius  and  in- 
versely porportional  to  the  speed  of  the  outer  flow.  At  600  miles  per  hour, 
this  would  correspond  to  t of  order  10“^  seconds  for  the  cylindrical  seg- 
ments characteristic  of  the  body  model  depicted  in  Figure  1.  Thus,  for  all 
intents  and  purposes,  it  may  be  postulated  that  separation  of  the  flow  takes 
place  at  all  limbs  almost  instantaneously  at  the  start  of  the  ejection  pro- 
cess near  sonic  speeds.  Moreover,  during  the  first  3 milliseconds  of  flow 
development,  the  separated  wake  stays  confined  to  the  immediate  region  of 
the  cylinder,  and  consists  essentially  of  two  thin  vortex-layers,  symmetri- 
cally situated.  It  is  thus  possible  to  model  this  situation  analytically  by 
simply  super-imposing  stationary  vortex  filaments  upon  the  linear  solution 
for  invlscid  flow  past  a cylinder  (see  Milne-Thomson,  1960,  pages  367-369 
and  Mello,  1959).  What  one  obtains  is  a reasonable  estimate  of  the  pressure 
distribution  along  the  surface  of  the  cylinder,  and  pressure  forces  are 
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those  of  greatest  concern  during  the  catapult  phase  of  lilgh-speed  ejections. 
Going  back  to  the  assumption,  then,  one  must  conclude  that  it  is  feasible  to 
begin  this  analysis  by  examining  and  super-imposing  potential  flow  solutions 
for  the  geometries  shown  in  Figures  1 and  2. 

The  third  major  assumption  whicli  shall  be  inherent  in  this  preliminary 
formulation  is  that  the  flow  is  two-dimensional.  That  is,  end  effects  and 
conditions  at  the  joints  of  the  body  model  will  be  temporarily  neglected  and 
the  outer  flow  as  it  approaches  the  limb  will  not  be  considered  to  be  a 
function  of  x,  the  coordinate  along  the  limb  centerline.  As  can  the  assump- 
tion that  follows,  this  one  may  be  subsequently  modified  by  introducing 
"correction  factors"  to  account  for  three-dimensional  effects  (see,  for 
example,  Goldstein,  1938,  page  439,  where  he  discusses  cylinders  of  finite 
span).  In  any  case,  the  consequences  of  assuming  a relatively  large  aspect 
ratio  (length  to  average  diameter)  for  the  limbs  should  be  to  predict  some- 
what higher  drag  forces  than  are  actually  manifest.  Effects  at  the  ends 
of  finite-span  elements  tend  to  alter  the  separated  flow  patterns  in  such  a 
way  that  profile  drag  is  actually  decreased. 

The  fourth,  and  final  major  assumption  which  will  be  made  is  that  the 
flow  is  incompressible.  Like  the  neglect  of  turbulence,  this  assumption  is 
clearly  arbitrary  at  this  point  and,  considering  that  the  ejection  Mach 
number  is  near  unity,  probably  not  a very  good  one.  Nevertheless,  it  may  be 
possible  to  include  the  effects  of  compressibility  through  the  use  of  cor- 
rection factors  based  on  the  results  of  wind  tunnel  experiments.  There  is 
some  data  available  which  presents  various  coefficients  of  drag,  lift,  etc. 
as  functions  of  the  mean  free-stream  Reynolds  number  and  Mach  number.  Per- 
haps such  data  can  be  incorporated  into  this  analysis  at  a later  time.  For 
the  present,  assuming  the  flow  to  be  incompressible  certainly  simplifies  the 
mathematics . 


SECTION  IV 

THE  COAXAL  COOPJIINATE  SYSTEM 


In  the  remainder  of  this  report,  we  shall  show  the  feasibility  of 
analyzing  the  aerodynamic  forces  exerted  on  an  articulated  body  subjected  to 
windblast  by  focusing  specific  attention  on  the  forearm  portion  of  the  right 
arm.  This  body  segment  was  discussed  in  some  detail  (Section  III)  with 
reference  to  Figures  1 and  2,  a.r.j  we  now  propose  to  examine  the  limb  in 
three  positions: 

1)  resting  and  pressing  against  an  arm  rest,  such  that  there  is  a flat 
surface  of  contact  between  forearm  and  seat  rest  (this  configuration  is 
depicted  in  cross-section  in  Figure  4); 

2)  resting  but  not  pressing  against  an  arm  rest,  such  that  the  contact 
between  forearm  and  seat  rest  is  just  a straight  line  element  of  both  sur- 
faces (this  configuration  is  Illustrated  in  Figure  1,  and  magnified  schemat- 
ically in  Figure  3);  and, 

3)  not  resting  at  all  against  any  surface,  such  that  the  forearm 
behaves  essentially  as  a slender  body  in  a free  stream  - disregarding,  for 
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the  moment,  its  attachment  to  the  other  parts  of  the  body  (this  configura- 
tion is  illustrated  in  Figure  2). 

The  above  three  positions  are  representative  of  some  which  are  typical 
during  ejection.  There  are  others,  but  they  have  not  as  yet  been  examined. 

In  order  to  pursue  the  analysis,  it  is  convenient,  at  this  point,  to 
introduce  a so-called  "coaxal"  coordinate  system  (see  Milne-Thomson,  1960, 
page  171),  depicted  in  Figure  4.  Consider  an  arbitrary  point,  P,  in  the  y-z 
plane.  The  location  of  P in  the  rectangular  coordinate  system  is  given  by 
the  point  (y,z),  while  in  the  cylindrical  coordinate  system,  it  is  the  point 
(r,9).  We  now  define  two  points,  G(-c,0)  and  E(+c,0)  on  the  y-axis  as 
shown,  and  locate  point  P as  the  intersection  of  the  lines  t2  and  r^^  drawn 

respectively  through  points  G and  E at  the  angles  62  and  0^.  Let  i = - 

0.,.  Then: 

sin  C = sin  (0^  - 62)  = sin  0^cos  ~ cos  0^sin  02 


. ^ cos  C y-c+2  . 2 2.2  ^ 

cot  5 = — : = T , and,  y - c + z - 2zc(cot  C)  = 0.  (14) 

sin  C 2zc 

2 2 

If  we  add  and  subtract  tlie  quantity  c cot  5 to  equation  (14),  it  may  be  put 
in  the  form: 

[z  - c(cot  Z) = c^(l  + cot^C)  = c^cosec^C*  (15) 

Equation  (15)  is  precisely  that  of  a circle  centered  at  the  point  y = 0, 
z = c(cot  C).  having  a radius  equal  to  c (cosec  C) • The  intersection  of  this 
circle  with  the  line  z = 0 is  the  chord  line  GE.  In  cross-section,  then. 
Figure  4 appears  to  be  a suitable  representation  for  case  1 described 
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earlier,  i.e.,  the  forearm  resting  and  pressing  against  an  arm  rest.  The 
forearm  cross-section  is  defined  by  the  circle  i = constant  and  the  flat 
surface  of  contact  between  forearm  and  seat  rest  (z  = 0)  is  defined  by  the 
chord  GE,  or  the  points  (+c,0),  (-c,0).  By  varying  c,  we  may  change  the 
extent  of  contact,  GE,  while  adjusting  ^ so  that  tlie  quantity  c (cosec  C) 
stays  within  physiologic  limits. 


Note  that  c(cosec  C)  = c(cot  C)  when  C = 0,  2ii , . . . 2mi  (n  = 0,  1,  2, 
...)  so  that  these  represent  the  limiting  case  of  the  circle  tangent  to  the 
y-axis  at  the  point  (0,0).  Furthermore,  when  C = 0,  c must  also  be  zero  to 
keep  c(cosec  C)  finite.  This  reduces  chord  GE  to  the  j^int  (0,0)  which  is 
consistent  with  the  fact  that  the  circle  becomes  tangent  to  the  y-axis. 

Thus,  C = 0 is  a convenient  way  to  represent  case  2,  i.e.,  the  arm  resting 
but  not  pressing  against  an  arm  rest.  Note  further  that  c(cot  C)  = 0 (for 

finite  c)  when  C = (2n  + l)-^  (n  = 0,  1,  2,  ...)  so  that  4 = represents 

the  limiting  case  of  the  circle  centered  at  the  point  (0,0)  of  radius  c.  We 
shall  show  later  that  this  is  a convenient  way  to  represent  case  3,  i.e., 
the  forearm  not  resting  against  any  surface  at  all. 


Let  us  now  define  a second  quantity,  n = Jin 


-"l/ 


Jin  r^  - Jin  r^,  so  that 


we  may  write:  n - = Jin  t2  + 102  - Jin  r^^  - 19^^  = Jin  (t2e  ) 


- Jin  (r^e  ),  or. 


Now,  r2e  may  be  written  as  (y  + c)  + iz  = (y  + Iz)  + c = x + c,  where  x 

101 

is  the  complex  variable  (y  + iz) . Similarly,  r^e  = r^cos  0^^ 

+ Tj^(i  sin  9j^)  = (y  - c)  + iz  = X “ c . Thus, 


n - iC 


The  coordinates  n«C  are  called  Coaxal  coordinates.  Equation  (15)  has  al- 
ready revealed  that  5 = constant  yields  a family  of  circles  centered  at  the 
point  y = 0,  z = c(cot  C)»  with  radius  c(cosec  5)»  c being  given.  It  can  be 
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shown  further  that  n = constant  yields  a family  of  circles  centered  at  the 

point  y = c(coth  n) , z = 0,  having  radii  equal  to  c(cosech  n).  Any  point  1' 

may  thus  be  defined  in  the  coaxal  coordinate  system  as  the  intersection  of 

two  circles,  C ‘ind  n • The  discussion  earlier,  and  the  material  which 
o o 

follows  indicate  that  the  use  of  coaxal  coordinates  is  convenient  for  tlie 
problem  being  considered  here. 

To  complete  this  section,  we  define  the  complex  variable  ^ = C + in. 
From  equation  (17),  then,  e (X  - c)  = x + c , from  which  we  get: 


= 

c 


e-i''’ 


+ 1 


e 


1 


(18) 


i 0 ^ 

Thus,  X ~ y iz  = re  = ic[cot  C = C + in  • From  these  relations,  one 

may  Invert  the  process  to  get  y and  z as  functions  of  c,  C and  ni  we  may 

—10 

also  define  the  complex  conjugate  functions:  x = y ~ iz  = re  = 

- icjcot  ^j,  C = C - in,  C = “ 02.  n = «.n 

■ SECTION  V 


THE  COMPLEX  VELOCITY  POTENTIAL 


For  the  two-dimensional  laminar  flow  of  an  incompressible  fluid,  the 
conservation  of  mass  requirement  may  be  identically  satisfied  by  introducing 
the  concept  of  a stream  function,  ip,  such  that,  in  a rectangular  cartesian 
coordinate  system, 

^ ■ I!  = If 

where : 

u = the  velocity  component  parallel  to  the  y-axis,  at  the  point  y,z, 

V = the  velocity  component  parallel  to  the  z-axis,  at  the  point  y,z, 

and  the  conservation  of  Mass  (Continuity)  equation  is  given  by: 


3u  ^ ^ 
9y  9z 


0. 


(20) 


If  the  fluid  is  considered  to  be  inviscid  as  well,  then  it  is  possible  to 
show  that  there  exists  a potential  function,  (j>,  such  that  the  velocity 

vector,  V,  has  components  corresponding  to  the  negative  gradient  of  4),  i.e.. 
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p 


u = - 1^  and,  V = - II  . (21) 

Substitution  of  the  expressions  (21)  into  the  Continuity  equation  (20) 
yields  the  well-known  Laplace  equation  for  the  determination  of  <p : 


2 2 

= 0 

2 2 
3y  3z 

Linder  the  assumptions  of  Section  III,  therefore,  we  may  introduce  a 
complex  velocity  potential,  w,  given  by: 


(22) 


and 


w = i()  + iij;  = w(x) 
w = - iij;  = w(x) 


(23) 


(2A) 


such  that: 


3w  dw  3x  dw  , . , 

- 3z  = - ^ 


dw 

dx 


= -u  + iv  (25) 


Thus , 


3w  _ dw  ^ dw  , . , 

+ — = + — + — (-i)  = -V  + lu, 

3z  3z 

dx  dx 


dw 

dx 


— = - u - iv  (26) 


dw  dw  2,2  2 

=u  + v = q 

dx  dx 


(27) 
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where  q is  the  square  of  the  magnitude  of  the  fluid  velocity  at  any  point. 
The  complex  velocity  potential  for  the  cross-flow  in  this  problem  may  be 
deduced  from  the  form  of  the  coaxal  coordinate  system  defined  by  equation 
(18)  of  Section  IV.  That  is,  consider  the  complex  potential: 

w = V cot  |=4i  + liJ;,  i;  = C + ih  (28) 

Expanding  this  potential,  and  separating  the  real  and  imaginary  parts 
yields,  after  some  algebra: 


(29) 
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2c 

n 2/£) 

1 , 2 /n\ 

tan  (— 

\nl 

' \n/ 

F 


Multiplying  the  top  and  bottom  of  equations  (29)  and  (30)  by  the  quantity 

2li\ 

cos  I— I reveals  that  <f)  remains  finite,  whereas  ij)  goes  to  zero  if  C = 0 (n 
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finite)  and  C = n^  (n  ^ 0).  Now,  the  limiting  stream  function,  ip  = 0 de- 
fines the  geometric  configuration  of  an  obstacle  disturbing  the  cross-flow 
of  velocity  V.  Thus,  consider  what  this  configuration  is  for  C = n^.  We 

have  already  seen  (from  equation  (15)  of  section  IV)  that,  for  c given,  C = 
constant  defines  circles  that  intercept  the  y-axis  at  the  points  +c . In 

particular,  then,  ^ defines  such  circles  for  various  values  of  the 

parameter  n (except  n = 0,  which  must  be  treated  as  a special  case  because  < 
and  p as  defined  by  equations  (29)  and  (30)  are  undefined  for  this  value  of 
n) . Furthermore,  the  degenerate  case,  C = 0,  n ^ 0 is  clearly  that  portion 
of  the  y-axis  lying  outside  the  points  +c , for  C = 9^  - 02  = 0 when  0^  = 6^= 

0 (positive  y-axis  to  the  right  of  point  +c  on  Figure  4)  and  C = 6^  ~ 02  “ 0 

when  0j  ~ ®2  ~ 1^0°  (negative  y-axls  to  the  left  of  point  -c). 


To  summarize,  we  find  that  the  complex  velocity  potential  (28)  defines 
exactly  the  cross-flow  over  a geometric  configuration  such  as  t) at  depicted 
in  Figure  4,  corresponding  to  case  (1)  of  the  forearm-seat-rest  interactions 
being  considered  here.  That  is,  the  limiting  streamline  if  = 0 defines  the 
arm-rest  y-axis  outside  the  points  +c,  together  with  the  circular  arm  inter- 
secting the  y-axis  at  the  points  +c.  We  may  point  out  further  that  the  same 
velocity  potential  would  also  define  the  flow  of  a fluid  past  a configura- 
tion which  consisted  of  that  shown  in  Figure  4,  together  with  its  mirror 
image  on  the  other  side  of  the  y-axis.  In  other  words,  if  a second  cylinder 
were  placed  on  the  negative  z side  of  Figure  4,  such  that  it  shared  a common 
base  GE  with  the  cylinder  shown,  then  the  complex  velocity  potential  (281 
would  characterize  that  flow  situation  as  well.  This  has  direct  relevance 
to  the  problem  of  one  limb  pressing  against  anotlier,  which  is  merely  an  ex- 
tension of  the  theory  being  developed  here. 


The  Special  Case  n = 0 

It  was  already  mentioned  in  Section  III  that  if  we  wish  to  keep  c(cosec 
5),  the  radius  of  the  circle  shown  in  Figure  4,  finite  as  C ^ 0,  then  c must 
also  go  to  zero,  reducing  the  chord  GE  to  the  point  (0,0)  and  making  the 
circle  taiigent  to  the  y-axis.  In  the  present  context,  this  means  that  c 

7T 

should  go  to  zero  as  n (and  hence,  C = trj)  does,  reducing  the  problem  to 

flow  around  a cylinder  tangent  to  the  y-axis,  which  corresponds  to  tlie 
cross-section  of  case  (2)  shown  in  Figure  3.  Note  from  equation  (18)  that 
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it'  X is  t<-’  remain  finite  as  c goes  to  zero,  tiien  tan  should  also  go  to 

zero  with  c,  indicating  that  C,  too,  becomes  very  small  with  n.  For  small 
c,  n and  we  may  therefore  write: 


= ic[cot  (I)]  = 


sin  ^ 


. 1 2ic 

1C  — = , or , C 

C C 


cfcosec  C)  = = — = — = a = radius  of  circle 

. IT  niT  niT 

sin  n-2  — j 


^10  7T  SL  ^ 

and  so,  = rrai,  - = = , which,  substituted  into  equation  (28), 

n n xn  X 

gives : 


w = V(Ttai)  cot  = iraV  coth  — = 4>  + iij^  , 

X X 


y + iz 


as  the  complex  velocity  potential  for  cross-flow  past  a cylinder  of  radius 
"a"  and  axis  parallel  to  x,  tangent  to  the  x-y  plane.  Once  again,  expanding 
equation  (33)  and  separating  the  real  and  imaginery  parts  yields  the 
velocity  potential  and  stream  function  for  the  motion: 


, airy  2 airz 
tanh  — ^ sec  — j 

„ r r 

,2  any  2 anz 

tanh  — 2 

r r 


anz  , 2 any 
tan  — 2 — 2 

r r 

'll  = naV. r q 

" , 2 any  2 anz 

tanh  — 2 — 2 

- r r 


where:  r'^  = y + z^  and  the  streamline  ij;  = 0 consists  of  the  boundary  z = 0 

2 2 

(the  real  axis)  and  the  circle  y + z = 2az  centered  at  z = a and  tangent 
to  the  real  axis  at  (0,0).  The  latter  is  obtained  in  analytic  form  from 
equation  (15),  using  the  approximations  (31)  and  (32). 
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Equation  ('35)  has  been  programmed  (Appendix)  and  Figure  5 shows  a 
typical  cross-flow  streamline  pattern  for  tlie  arbitrary  case  of  a = V = 1 ' 

whicli  essentially  non-dimens  ionalizes  the  stream  function  ip . Note,  in  ' 

particular,  the  stagnation  point  at  y =0,  z = 0 and  the  relatively  high  ' 

velocity  (proportional  to  the  distance  between  streamlines)  at  y = 0,  | 

z = 2.  This  indicates,  as  will  be  shown  more  rigorously  later  on,  that  the 
pressure  at  the  base  of  the  cylinder  is  much  higher  than  it  is  on  top  of  the 
cylinder,  so  that  there  is  a tendency  to  "lift"  it  from  the  surface  with 
which  it  is  in  contact.  Note,  also,  tliat  the  streamline  pattern  evens  out 
to  a nearly  parallel,  uniform  distribution  as  we  move  further  and  further 
away  from  the  cylinder.  In  fact,  if  we  take  the  limit  of  equation  (33)  as 
X we  find  that: 


TraV 


Vx 


(36) 


X 


whicli  is  exactly  the  complex  velocity  potential  for  a uniform  stream,  moving 
from  right  to  left  parallel  to  the  real  axis.  As  an  interesting  side-light, 
we  might  also  mention  that  the  complex  potential  (36)  can  be  deduced  from 
equations  (18)  and  (28)  as  the  special  case  for  n = 2,  so  that  the  latter 
represents  infinitely  unbounded  flow  parallel  to  the  real  plane  in  the 
absence  of  an  obstacle.  This  has  no  immediate  relevance  to  the  problem  at 
hand,  but  constitutes  an  interesting  observation. 


The  Special  Case  n = 1 

For  n = 1,  the  complex  potential  (28)  reduces  to: 

w = 2ciV  cot  C = 2ciV  ^ . (37) 

sin  c 


Making  use  of  the  trigonometric  identities  for  the  sine  and  cosine  of  double 
angles,  we  may  re-write  equation  (37)  as  follows: 


w = 2ciV 


2 5 . 2 r; 

cos  y - sin  2 

2sin  -2  cos  2 


= v|^ic  ^cot  - Ic  ^ 


tan  j 


ic^cot  y) 


ic  ^cot 


(38) 
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0-255 


(39) 


\ 


But 


, from  equation  (18),  ic^cot  = X>  so  equation  (38)  reduces  to: 


w = V 


'1'  + iij'  . X = y + iz  = re 


which  is  precisely  the  complex  velocity  potential  for  uniform  streaming  past 
a free  cylinder  of  radius  c oriented  with  its  axis  perpendicular  to  the 
flow.  This,  then,  corresponds  to  the  cross-flow  at  a station  x of  case  (1) 
which  we  are  examining  and  which  is  depicted  in  i'igure  2 (note  that  the 
velocity  V is  perpendicular  to  the  cone  axis,  x) . It  is  now  clear  that  the 
coaxal  coordinate  system,  and  the  complex  potential  (28)  are  particularly 
useful  for  describing  the  cross-flow  in  each  of  the  three  configurations 
under  investigation  in  this  study.  In  fact,  they  can  be  used  to  describe  a 
variety  of  other  configurations  as  well  - for  example  - one  limb  in  contact 
with  another  as  discussed  at  the  end  of  the  section  dealing  with  the  case 
n - 0. 


As  has  been 
separate  it  into 


done  in  previous  sections,  we  now  expand  equation  (39)  and 
its  real  and  imaginery  parts  to  obtain: 


4,  = 


<P  = 


e 


0 


(40) 


(41) 


At  tliis  point,  the  complex  velocity  potential,  stream  function  and  potential 
function  have  been  deduced  for  each  of  the  three  basic  configurations  de- 
picted In  Figures  2,  3 and  4.  They  are,  respectively,  equations  (28),  (29) 
and  (30)  for  Figure  4,  corresponding  to  case  1 (forearm  pressing  against 
seat  rest  handle);  equations  (33),  (34)  and  (35)  for  Figure  3,  corresponding 
to  case  2 (forearm  just  touching,  i.e.,  tangent  to  seat  rest  handle);  and, 
equations  (39),  (40)  and  (41)  for  Figure  2,  corresponding  to  case  3 (forearm 
freely  exposed  to  oncoming  stream).  In  the  section  that  follows,  these 
velocity  functions  shall  be  used  to  estimate  the  pressure  distribution 
around  the  forearm  in  each  configuration,  and  we  shall  see  how  the  angle  of 
attack,  a,  and  the  cone  half-angle,  8,  enter  into  the  analysis. 


SECTION  VI 

ESTIMATION  OF  THE  PRESSURE  DISTRIBUTION  AROUND  THE  FOREARM 


Ihe  stream  function,  , and  velocity  potential,  (f , liave  both  been  shown 
to  satisfy  the  conservation  of  mass  equation  (20).  In  order  to  uniquely  de- 
fine the  fluid  motion,  hciwevt i , we  must  also  satisfy  conservation 
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relationships  for  momentum  and  energy.  The  latter  two  become  uncoupled  from 
one  another  under  the  assumption  of  incompressibility,  and  it  thus  becomes 
possible,  for  an  inviscid  fluid,  to  solve  directly  for  the  pressure  distri- 
bution once  <p  is  obtained  from  equation  (22).  Tlie  procedure  involves  the 
solution  of  an  integrated  form  of  the  conservation  of  momentum  equation, 
where  the  integration  is  performed  along  a streamline  (ip  = constant)  of  the 
flow.  T(ie  result  is  the  well-known  Bernoulli  equation,  which,  neglecting 
body  forces,  takes  the  form: 


,,  2 

- ^ ^ + ^ = Constant,  C (A2) 

9t  2 p o 

where ; 

())  = the  local  cross-flow  velocity  potential,  determined  from  equation 

(22); 

t = time; 

q^  = the  magnitude  squared  of  the  local  fluid  velocity  at  a specified 
point  in  the  flow,  as  determined  by  equation  (27)  utilizing  the 
complex  velocity  potential; 

p = the  pressure  due  to  cross-flow  at  the  same  point  as  that  at  which 
q^  is  determined,  for  a corresponding  time,  t; 

p = the  fluid  mass  density;  and, 

C = a constant,  generally  determined  from  known  ambient  conditions, 
or  from  conditions  at  a known  point  in  the  flow. 


The  integrated  momentum  equation  (42)  is  obviously  a form  of  energy  equa- 


pressure  energy  term,  p (proportional  to  molecular  kinetic  energy  per  unit 
mass  of  fluid) , and  so  on  (potential  energy  terras  also  appear  if  body  and 
field  forces  are  taken  into  account  in  the  momentum  equation).  Thus,  the 
constant  represents  the  "total"  energy  of  the  flow  along  a streamline, 

which  remains  constant  in  the  absence  of  dissipating  viscous  effects. 


Equation  (42)  may  be  used  to  estimate  the  pressure  distribution  for 
each  of  the  three  configurations  shown  in  Figures  2,  3 and  4.  Thus,  con- 
sider the  complex  velocity  potential  given  by  equation  (28)  , together  with 

— 2ci  c 

its  complex  conjugate,  w = -V cot  — . From  equation  (27)  we  have: 


2 

q 


dw  dw 


dw  d^  dw  dg 


2 

16V^ 

4 

n 


I sin 
sin 


C 7 2 

2 2 I 

^sin^l 


(43) 
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wliere  the  derivatives  are  obtained  from  equations  (28),  (18)  and  their  res- 
pective conjugates.  Expanding,  rearranging  and  simplifying  equation  (43) 
yields : 


2 _ 16V^  I cosh  n ~ cosh  C 


q = 


cosh^n 


(44) 


- cos 


which  is  the  second  term  of  equation  (42).  We  can  also  use  the  chain  rule 
of  calculus  to  write: 


8(|)  _ ^ ^ ^ d2^ 

3t  9c  9t  9c  dx  dt 


From  equation  (29) : 


(45) 


) 

|o? 

1 2c 

'n' 

Isec^l 

fi' 

1 1 

9c  9c  j 

1 n 2, 

1 tan  1 

(nj 

I + tanh^l 

(46) 


Also,  by  definition: 
dx 

— = - U = -U^  cos  a = the  stream  velocity  parallel  to  the  x-axis. 


dc  dn  , dc  , 


dc 

The  quantity  — defines  the  surface  of  contact  between  the  forearm  and  rhe 

seat  rest.  It  is  finite  when  the  forearm  configuration  is  taken  to  be  a 
truncated  cone,  and  its  value  depends  upon  6,  the  cone  half-angle,  and  n, 
which  determines  the  "degree"  or  extent  of  pressing  (i.e.,  contact)  between 

d c 

forearm  and  arm  rest.  In  fact,  we  may  write: 

dx 

indicates  how  the  chord  line  GE  (Figure  4)  varies  with  n - hence  the  extent 

of  contact  - and  ^ = f(0)  Indicates  how  the  chord  line  GE  varies  with  x - 

hence  the  half-angle  of  the  forearm  truncated  cone  configuration.  Specifi- 
cation of  these  quantities,  together  with  U^,  the  angle  of  attack,  a,  and 
equation  (46)  allow  equation  (45)  to  be  evaluated,  which  then  determines  the 
first  term  of  equation  (42).  Adding  appropriate  ambient  information  allows 
to  be  computed.  Thus,  combining  equations  (44),  (45)  and  (42)  with 

known,  we  may  solve  for  p,  the  pressure  distribution  on  the  surface  of  the 
forearm  as  a function  of  several  variables  - which  include  forearm  geometry, 
degree  of  forearm/seat  rest  contact,  free-stream  conditions  and  the  manner 
in  which  the  stream  approaches  the  limb.  Obviously,  the  solution  process 
becomes  quickly  complicated  and  computer  methods  must  be  employed  to  handle 
the  mathematics  in  a convenient  fashion.  For  illustrative  purposes. 
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Iiowever,  we  may  go  directly  to  two  specific  cases  in  order  to  show  the 
feasibility  of  the  approach  and  continue  the  analytic  development. 


The  Special  Case  n = 1 


The  easiest  case  to  examine  first  is  case  3 (n  = 1)  of  Figure  2,  where 
the  forearm  is  depicted  as  a free,  tapered  slender  body  at  an  angle  of 
attack  to  an  oncoming  stream.  Thus,  consider  the  complex  velocity  potential 

- f- 

given  by  equation  (39),  together  with  its  complex  conjugate,  w = V x + ~ • 
From  equation  (27),  we  have:  ^ 


2 dw  dw  ,,  , c 

q = q =V1 ;r 

dx  2 

dx  X 


1 - — 2 , or,  expanding  and  rearranging 


2 „2 
q = V 


\(  2 / c2\^  2I 

( 1 - 1 cos  0 + (1  + -2  ) sin  6 . 


Equation  (45)  now  becomes,  .since  c = R(x)  in  Figure  2: 


3tj)  _ ^ ^ ^ ^ ^ 

3t  3c  3t  3R  dx  dt 


From  equation  (40): 


I!  ■ '[t] 


cos  6 = , since  c = R in  case  3. 


Also,  from  Figure  2: 


UR  „ / . . . 

— - = tan  g (a  negative  quantity). 


And,  by  definition: 


— = -U=-U  cos  a = the  stream  velocity  parallel  to  the  x-axis. 
dt  o 


i'uit  ing  all  of  this  information  into  equation  (48),  we  get; 
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= U cos  0 tan  3 sin  2a  (49)  | 

Note  in  equation  (49)  that  two  negative  quantities  exist,  but  tiieir  sign  is  | 

inherent,  so  that  we  need  concern  ourselves  only  with  their  absolute  values.  j 

That  is,  the  complex  potential  <}>  defined  by  equation  (40)  is  for  a stream  j 

moving  with  velocity  V in  the  negative-y  direction,  and  the  fact  that  R(x)  j 

decreases  with  x has  been  accounted  for  in  the  first  step  above  by  putting  a j 

minus  sign  in  front  of  tan  3-  Should  the  circumstances  of  Figure  2 ciiange 
from  those  indicated,  care  must  be  taken  to  make  the  appropriate  sign 
changes  where  necessary.  Essentially,  equation  (49)  postulates  an  analogy 
between  the  variation  of  cross-flow  with  x on  a truncated  cone  configura- 
tion, and  the  variation  of  cross-flow  with  t on  a circular  cylinder  config- 
uration started  impulsively  from  rest  in  a stationary,  unbounded  fluid  (for 
further  discussion  of  this  point,  see  Allen  and  Perkins,  1951). 

Now,  if  we  let  V = U sin  a in  equation  (47)  and  substitute  (47)  and 

° j 

(49)  into  (42),  we  get:  j 


- U 


^(f) 


cos  6 tan  3 sin  2a  + -^  = C 

P o 


(50) 


As  r ->■  “>,  p approaches  some  free  stream  pressure,  p^,  so  that: 


^ ‘^0^1,2  .2 
C = — + — U sin  a 
o p 2 o 


(51) 


Equation  (50),  with  defined  by  equation  (51)  can  be  used  to  calculate  the 

pressure,  p anywhere  in  the  flow  field.  We  are  most  Interested,  however,  in 
what  happens  on  the  surface  of  the  cone.  Thus,  let  r = R to  get: 


V ^ sin'^  a (4sin^  6)  - U cos  9 tan  3 sin  2a  + -^  = — + 4 1.'  ^ sin^a  (52) 

/ C\  C\  / 


2 o 
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sure 


Ac  this  point, 
coefficient,  C 


it  is  convenient 

p - p 

Then 


o 


1 n 2 
0 pL' 

2 o 


to  define  a dimensionless 
, from  equation  (52); 


surface  pres- 


C 

P 


. 2 

sin  a 


1-4 


sin^  0| 


+ 2 cos  0 tan  3 sin  2a 


(53) 


Note,  in  equation  (53)  that  C =0  (i.e.,  p = p ) if  a = 0°,  that  is,  if  the 

p o 

stream  approaches  the  limb  parallel  to  its  centerline,  the  x-axis.  In  fact, 
this  is  strictly  true  only  if  3 = 0°  as  well,  i.e.,  if  the  limb  configura- 
tion is  a perfect  right  circular  cylinder.  If  3 0,  such  that  the  limb 

configuration  is  that  of  a truncated  cone  facing  a stream  coming  along 
parallel  to  the  x-axis,  then  the  surface  pressure  at  any  station,  desig- 
nated p _,  will  differ  slightly  from  the  static  pressure,  p , for  small  6- 

The  discrepancy  in  equation  (53)  stems  from  the  fact  that  the  derivation  of 
this  equation  for  C has  concerned  itself  only  with  the  effects  due  to  a 

cross-flow,  V,  coming  at  the  limb.  Account  has  not  been  taken  of  the  pres- 
sure distribution  generated  on  the  inclined  conical  surface  as  a result  of 
the  axial-flow  component,  U.  In  the  work  of  Allen  and  Perkins  (1951),  an 
argument  is  presented  for  considering  tliat  this  latter  pressure  distribu- 
tion remains  constant  with  a,  adds  linearly  to  C for  small  3 and  reduces  to 

P 

zero  for  3=0.  Thus,  we  may  simply  add  an  arbitrary  correction  factor  to 
equation  (53),  in  the  form  of  a pressure  coefficient  designating  the  effects 
on  the  surface  due  to  U: 


C = sin^  a [1 
P 


4 


. 2 
sin 


6] 


+ 2 cos  0 tan  3 sin  2a  + C 

Pa=0 


(54) 


where  it  is  assumed  that  C does  not  vary  much  with  a (for  small  3)  and 

Pa=0 

goes  to  zero  with  3.  For  the  time  being,  this  latter  pressure  coefficient 
needs  to  be  determined  empirically,  and  we  shall  neglect  it  temporarily.  In 
fact,  it  turns  out  that  the  contribution  of  C to  total  cross-force  cal- 

Pa=0 

culations  is  not  important  because  its  integral  around  the  surface  vanishes 
in  this  direction  (Allen  and  Perkins,  1951). 


Note  further,  that,  when  3=0 

2 

(sin  a = 1,  sin  2a  = 0)equation  (54) 
the  pressure  distribution  on  a right 


[C  = 0,  tan  3=0]  and  a = 90° 

^a=0 

reduces  to  the  classical  equation  for 
circular  cylinder  oriented  with  its 


axis  normal  to  a stream  coming  at  it  with  velocity  V 


(tan  a = - = 


U = 0) 
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C =1-4  sin^  e (55) 

P 

Kor  tlie  problem  under  consideration,  large  values  of  a are  not  unlikely  be- 
cause the  presence  of  a torso  (Figure  1)  defiecting  an  oncoming  stream  of 
fluid  would  cause  that  fluid  to  approach  tlie  forearm  at  a ratlier  large  angle 
of  attack. 


Equation  (53)  is  shown  plotted  in  Figures  6 and  7 for  two  values  of  5, 
with  a acting  as  a parameter.  These  two  figures,  together  wi  th  equations 
(47),  (52)  and  (53)  may  be  used  to  illustrate  several  important  points. 
First,  we  observe  in  equation  (47)  that  stagnation  points  (q  = 0)  exist  on 
the  surface  (r  = c = R)  of  the  limb  at  9 = 0°  and  180°.  The  stagnation 
point  at  0°  is  undoubtedly  correct,  whereas  the  one  at  180°  must  be  viewed 
more  carefully  in  terms  of  the  separated  flow  patterns  on  the  lee  side  of 
bluff  bodies.  For  9=0°,  equation  (53)  becomes: 

2 

C = sin  a + 2 tan  8 sin  2a  (56) 

P 

1 2 

and  we  see  that  C =1  for  a - 90°  (p  = p + — pU  = maximum) , C =0  for 
p o 2 o p 

a = 0°  (p  = Pq) • Thus,  0 ^ Cp  £ 1 for  0°  £ a 90°  with  9=0°  (or  180°), 
which  is  shown  clearly  in  Figures  6 and  7. 


Second,  if  we  differentiate  equation  (52)  with  respect  to  9 and  set  the 
result  equal  to  zero,  this  yields: 
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from  which: 


^ = p|  ^_^sin^a(8sin  9 cos  9)  (-1)  + (-sin  9)  tan  6 sin  2a  | = 0 (57) 


aj  = 0 


cos  9 = - 


tan  8 1^ 
tan  a 2 


(58) 


Observe,  then,  that,  for  8 finite  (not  equal  to  zero),  cos  9 varies  as  the 
inverse  of  tan  a.  Substituting  equation  (58)  into  (53)  gives: 


tan  8 sin  2a 


, 1 


.2  I,  , ^ 2 „ cos 
= sin  a|l  + tan  8 ■ - y 

sin 


_ 4'  - 

a 


2 ^ . ,cos  a 

tan  6 (2sin  a cos  a) — ;; 

sin  a 
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tan^cose  sLn  ( l~4sLn^e 


Figure  7.  Surface  pressure  distribution  around  lore.irm  having  four  degree  angle  of  taper 
and  sitting  free  in  wind  stream. 
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It  which  the  pressure,  p, 
and  the  cone  half-angle, 


has  its  minimum  value  based 

point  C = -3  (c.f.,  equa- 
H P 

t ion  (59)  and  Figures  6 and  7).  In  fact,  we  see  (from  equation  (59))  that 
this  is  true  for  a = 90°  regardless  of  the  value  of  3,  and  equation  (58) 
shows  tliat  in  this  case  cos  6 = 0,  so  that  0 = 90°.  Now,  for  the  angle  of 
attack  decreasing  from  90°,  the  value  of  tan  a decreases,  causing  cos  0 to 
become  more  negative  in  equation  (58)  and  thus  forcing  0 to  increase  from 
90°.  Physically,  this  means  that  the  point  of  minimum  pressure  moves  more 
lee-ward  (i.e.,  towards  the  rear  of  the  cone  facing  the  flow)  as  the  angle 
of  attack  decreases.  Since  the  flow  tends  to  separate  off  of  a bluff  body 
near  the  point  of  minimum  pressure,  it  follows  that  the  point  of  flow 
separation  moves  lee-ward  with  decreasing  a for  the  geometric  configuration 
depicted  in  Figure  2.  Such  motion  depends  weakly  on  the  angle  6 (for  3 
small)  in  accordance  with  equations  (58)  and  (59)  above,  and  it  is  illus- 
trated in  Figures  6 and  7.  Here  the  minimum  value  of  the  curves  shown  is 
seen  to  be  drifting  slowly  to  the  right  (increasing  0)  as  ci  decreases  from 
90°.  Note  also  in  these  figures  that  the  pressure  minimum  increases  with 
decreasing  angles  of  attack,  such  that  there  is  less  pressure  to  be 
recovered  on  the  back  side  of  the  limb  as  its  axis  comes  more  in  line  with 
the  incident  flow.  In  the  limit  of  a -*■  0 , there  is  no  pressure  loss  at  all 
as  the  fluid  traverses  0.  Tliat  is,  C = 0 as  previously  discussed  with  re- 
spect to  equation  (53).  However,  recall  that  this  is  strictly  true  only  if 
we  can  neglect  C , which  we  have  done  temporarily  in  the  present  invest i- 
^a=0 

gation.  Finally,  observe  that  the  curves  shown  in  Figures  6 and  7 for  a = 
90°  are  exactly  those  of  equation  (55),  which  is  the  standard  right-circu- 
lar-cylinder solution. 


To  briefly  summarize,  then,  the  analytic  results  of  this  section 
reveal : 


1)  that  pressure  maxima  occur  at  0=0°,  corresponding  to  points  of  flow 

stagnation  for  a ^ 0°  --  and  that  these  maxima  lie  between  p and  p 

o *^0 

1 2 

+ Y pU^  , depending  on  the  angle  of  attack; 

2)  that  pressure  minima  occur  between  0 = 90°  and  0 = 180°,  depending  on  a, 
and  correspond  roughly  to  points  near  where  the  flow  will  separate  Tor 

o 3 2 

Cl  / 0°.  The  pressure  at  these  points  lies  between  p and  p - -r-  pF  , 

o o 2 o 
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A 


and  Lilt'  points  LliomsoJ  vos  move  It-oward  ( i no  rt-as  i ng  Oj  with  decreasing  a 
at  a rate  mildly  dependent  on  3 (lor  3 small). 

The  Special  Case  n = 0 

We  move  on,  now,  to  examine  the  case  of  the  forearm  resting,  but  not 
pressing  against  an  arm  rest  (case  2,  Figure  J,  n = 0).  Thus,  consider  the 
complex  velocity  potential  given  by  equation  (3'3)  , together  with  its  complex 

conjugate,  w = TtaV  coth  . From  equation  (27),  we  have; 

X 


2 dw  dw  a IT  V , 2 air  , 2 an 

q = -r-  — = — ■-  cosech  — cosech  — 

dx  2—2  X ~ 

dX  X X X 


(bO) 


which,  upon  expanding  in  accordance  with  the  definition  of  x>  X>  and  per- 
forming some  algebraic  simplifications,  becomes: 

I 

2 

, 2airv  2anz 

cosh  — ^ - cos  — — 

r r 

2 2 2 2 
where:  r = y + z . On  the  surface  of  the  limb,  r = 2az  (from  equations 

(15)  and  (32)),  so  that  equation  (61)  becomes: 


2 a IT  V 

q = ^ 


2 a IT  V , 4 TTv 
q = — -y-  sech 

4z 


(62) 


Similarly,  on  the  plane  to  which  the  cone  is  tangent,  z = 0,  so  that  equa- 
tion (61)  reduces  to: 


2 a IT  V , 4 an 

q = : cosech  — . 

_ q y 


For  the  case  n = 0,  equation  (45)  becomes. 


3t 


9 a 9t 


9<f  ^1  ^ 
9a  dx  dt 


(63) 


(64) 


where  is-  given  by  equation  (34)  and  "a"  is  the  radius  of  the  forearm 
(equation  (32)),  which  varies  from  the  wrist  to  the  elbow  (along  x)  as 
shown  in  Figure  3.  From  equation  (34): 
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2aTTz 
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aTTz\  airy  ( , 2 atiy  2 airz\ 
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(65) 


Again,  we  may  write; 


and , 


da  „ / . . N 

— = tan  p (a  negative  quantity), 


(66) 


dx 

— = - U cos  a . 
dt  o 


(67) 


To  simplify  the  mathematics,  recall  tliat  we  are  most  interested  in  what 

2 

happens  on  the  surface  of  the  limb.  Thus,  substituting  r = 2az  into  equa- 
tion (65)  and  putting  that,  together  with  the  expressions  (66)  and  (67)  into 
(64) , yields . 


3i  = 
at 


tt(V  si 
o 


n ot)  1^1 


2 IZ1 


tanh  ^ ^ sech^  ^1  (-tan  6)  (-U  cos  a)  (68) 

2z  2z  2zJ  o 


We  comment  here  that  equation  (65)  needs  to  be  multiplied  by  cos  — — on  top 

r 


and  bottom  in  order  that  its  limit  as  r 


2az  may  be  determined.  Further- 
da 


more,  a negative  sign  has  again  been  inserted  before  the  ^ term  in  equation 

(68)  in  order  to  account  for  the  fact  that  tan  S is  a negative  quantity  as 
discussed  previously.  Finally,  it  is  a relatively  straightforward  matter  to 
show  that  the  limit  of  equation  (65)  as  r is  zero,  and  that  of  equation 

2 

(61)  is  V . Thus,  from  Bernoulli's  equation  applied  at  infinity,  we  again 
arrive  at  the  conclusion  that; 


r - “ ^ 2 ,T  2 . 2 

C = — + — U sin  a 
op  2 o 


(69) 


and  therefore,  equation  (42),  written  for  the  surface  of  the  forearm  using 
equations  (62),  (68)  and  (69),  becomes; 
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^ = 1-  ■;r  U sin  a 

P P 2 o 


",  in  the  form  of  the  non-dimensional ized  surface  pressure  coefficient. 


C = ; = sin^  a 1 - — ^ sech"*  + usin  2oi  tan  B ftanh  ^ 

P |pu/  L 4z^  2zJ  L 


+ ^ sech 
2z 


2 7^1 

2zJ 


«ihere  ^ = cot  0 (see  Figure  A)  , r = 2az  and  we  are  still  neglecting  the 
z 

contribution  from  the  axial  flow  component,  U,  in  the  form  of  C as  we  did 

P(i=0 

before.  That  is,  note  again  in  equation  (71)  that  = 0,  i.e.,  p = when 
a = 0°.  Now,  from  Figure  4,  it  is  possible  to  write; 

2 2 

sin^  = I — = ■—  , so  that:  z^  = 4a^sin^  6 (72) 

2 2az  2a 
r 

Substituting  equation  (72)  into  equation  (71),  we  may  write: 


Ti  sech  0) 

C = sin^  a 1 q; + usin  2atan  B 

P 16sin  0 


j^tanh^-^ot  ej 


,11  „ ,2 

+ ^ot  0 sech 


(^ot  0) 


which  is  an  expression  annlagous  to  equation  (53),  giving  as  a function 
of  a,  B,  and  0. 

Wlion  6 = 0°  fc  =0,  tan  P = o]  and  a = 90°  (sin^  n = 1,  sin  2a  = 0) 


equation  (73)  reduces  to: 


(74) 


C 
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4 4/  Ti 

TT  seen  ly  cot 
_ 

16sin  6 


whicli  is  the  equivalent  of  equation  (55)  written  for  a free  rigtit  circular 
cylinder.  It  is  interesting  to  compare  the  results  obtained  from  equations 

(55)  and  (74).  For  example,  at  the  top  of  the  cylinder  ^9  = equation 


(55)  gives  Cp 


-3,  whereas  equation  (74)  yields  C 


-5.  Thus,  the  pres- 


sure coefficient  in  the  latter  case  is  some  67%  lower  than  it  is  for  the 
free  cylinder,  indicating  that  tlie  presence  of  the  flat  surface  causes  the 
pressure  at  the  top  of  the  cylinder  to  be  greatly  reduced  compared  with  the 
corresponding  case  where  the  surface  of  contact  is  absent.  Moreover,  for  6 

2tt 


small,  cot  0 = "T 

V 


sech 


16e 


and  sin  6 = 0 , so  that  equation  (74) 


becomes : 


C = 
P 


(75) 


the  limit  of  which,  as  6 ->  0 is  1 . At  the  bottom  of  the  cylinder,  where  it 
is  in  contact  with  the  flat  surface,  we  therefore  have  a stagnation  point 
(q  = 0)  at  which  = 1 for  a = 90°.  The  pressure  at  this  point  is  much 

higher  than  it  is  at  the  bottom  of  a cylinder  sitting  free  in  an  unbounded 

flow,  where  equation  (55)  reveals  that  = -3  when  ^ ~ ~ Y'  effect  of 

the  flat  surface,  then,  is  to  distort  the  symmetric  pressure  distribution 
around  the  cylinder,  such  that  the  pressure  is  higher  at  the  bottom  (9  = 0° 
in  Figure  4)  than  it  is  at  the  top  (0  = 90°).  In  fact,  it  is  higher  by  an 
amount : 


0=0 


^S^e=90 
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32 


(76) 


There  is  thus  a tendency  for  the  cylinder  to  be 
maximum  pressure  gradient  of  some  six  times  the 

1-  I,  2 

sure,  -T  pU 
2 o 


lifted  from  the  surface  by  a 
free  stream  dynamic  pres- 
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At  the  stagnation  point  0 = 0,  where  the  static  pressure  lias  its  maxi- 
mum value  on  the  surface  of  the  conical  limb,  equation  (73)  reduces  to: 


C = sin"  a + Ti  sin  2a  tan  6 (77) 

P 


which  is  similar  to  equation  (56),  differing  only  in  the  coefficient  of  the 
second  term  on  the  right  hand  side.  The  first  term  follows  directly  from 
the  discussion  which  has  preceded  equation  (75).  The  second  teem  results 

from  the  fact  that,  for  9 small,  tanh^-^otjs  + ^ot  6 sech^  ~ 


TT 


tanh 

Thus 


, again  we  see  that  C = 


IT 

26 

1 


for  a 


the  limit  of  which  is  1 as  9 0. 

= 90°  f'  = p + 4 pU  ^ = maximum),  C = 
\ o 2 o ' P 


0 for  a = 0°  (p  = p ),  and  0 < C <1  for  0°  < a < 90°.  The  difference  be- 
"^o  — p — — — 

tween  this  case  and  the  previous  section,  however,  is  that  equation  (77) 
holds  at  the  bottom  of  the  forearm,  whereas  equation  (56)  holds  at  the  for- 
ward stagnation  point  along  the  front  side  of  the  limb. 


At  the  top  of  the  cone  equation  (73)  becomes  (for  0 = 90°) : 


2 

C = sin  a 1 

P L 


(78) 


so  that: 


P,  - P 6 

\ /n  \ _ bottom  top  . T ^ 

^^p^e=0  ~ ^^p^0=9O  12  ^ 16 

2 f^^o 

The  pressure  gradient  tending  to  "lift"  the  forearm  off  of  the  seat  arm  rest 
thus  appears  to  increase  with  a up  to  a maximum  for  a = 45°.  However,  the 
variation  of  Ap  with  a is  scaled  by  tan  S,  which  is  very  small  in  the  physi- 
ologic range.  For  6 as  high  as  10°,  for  example,  the  maximum  value  of 

4 

IT  sin  2a  tan  6 is  only  0.55,  compared  with  the  value  of  — which  is  6.09. 

16 

This  may,  or  may  not  be  significant  for  very  large  values  of  U^,  and  viewed 

in  terms  of  the  ability  of  the  forearm  musculature  to  resist  displacing 
forces,  it  could  turn  out  to  be  critical.  In  any  event,  the  tendency  for 
limb-displacement  seems  to  be  substantial  in  cases  where  the  interaction  of 
the  limb  with  a solid  surface  produces  a stagnation  point  in  the  flow  — 
especially  if  the  flow  is  high-speed. 

If  we  differentiate  equation  (73)  with  respect  to  9 and  set  the  result 
equal  to  zero,  this  is  equivalent  to  finding  the  point  along  the  surface 
where  the  pressure  as  defined  by  equation  (7T)  is  a minimum.  Letting 

6 = Y cot  9,  it  is  found,  after  some  algebra,  that  9 must  satisfy  the 
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re  Lat ion : 


2 

sin  20  - ir  tanh  6 secli  6 _ ^ tan  ji 
sin  20  tanh  6+2  . 4 , 2 tan  u 

sin  6 71 


(80) 


wliicti  is  .inalajtous  to  equation  (58)  for  the  case  of  the  free  cone.  Thus, 
again  we  find  that  the  angle  for  wliich  the  pressure  is  minimum  varies  with 
both  a and  0,  wltli  a = 90°  causing  to  have  its  minimum  value  (just  under 

-5)  at  0 = 90°  — the  very  top  of  the  limb.  This  is  clearly  illustrated  in 
Figures  8 and  9,  wliere  equation  (73)  is  plotted  for  two  values  of  6,  with  u 
acting  as  a parameter.  Comparison  of  these  two  Figures  with  Figures  6 and  7 
reveals  a qu.ilitative  similarity  in  the  ,'eneral  character  of  the  curves 
shown,  but  aiso  highligh*j|  some  dramatic  differences.  For  example,  note  how 
flat  the  curves  ot  Figures  8 and  9 are  for  0 between  0°  and  40°,  and  again 
between  140°  and  180°.  This  represents  a very  high  pressure  region  (C^  > 0) 

spanning  some  160°  ot  arc  along  the  lower  half  of  the  circular  cross  section 
in  contact  with  the  flat  surface  (note  that  an  inscribed  angle  has  the  same 
measure,  in  angle  degrees,  as  half  of  its  intercepted  arc,  in  arc  degrees). 

For  6 between  40°  and  140°  the  pressure  drops  off  quite  sharply,  and  there 
results  a significant  low  pressure  region  for  the  remaining  200°  of  arc. 

Thus,  Figures  8 and  9 reinforce  the  concept  that  large  forces  exist  which 
tend  to  dislodge  a limb  from  a surface  with  which  it  is  in  contact.  In  the 
present  context,  such  forces  arise  as  a result  of  causing  a high-speed  flow 
to  stagnate.  Later  on,  we  shall  discuss  (in  principle,  for  the  moment) 
additional  forces  that  arise  as  a result  of  causing  a flowing  stream  of 

fluid  to  separate  from  (or  to  fail  to  follow  the  contour  of)  a solid  surface  j 

with  which  it  is  in  contact.  Note  further  in  Figures  8 and  9 that  the  i 

minima  of  the  curves  shown  also  drift  to  the  right  (increasing  9)  with  de- 
creasing angle  of  attack,  and  that  their  magnitude  increases.  As  was  the  i 

case  for  the  free  cone,  this  means  that  the  point  of  minimum  pressure  moves  j 

lee-ward  as  the  angle  of  attack  decreases  (as  does  the  corresponding  point  j 

of  flow  separation)  and  there  is  also  less  pressure  to  be  recovered  on  the 
back  side  of  the  limb  as  its  axis  comes  more  in  line  with  the  incident  flow. 

These  events  are  dependent  on  the  angle  0 in  accordance  with  equation  (80). 

Summarizing,  then,  the  essential  analytic  findings  of  this  section,  we 

see : 

1)  that  pressure  maxima  occur  at  the  point  of  contact  between  a conical 

surface  and  a flat  restraining  plate,  varying  in  magnitude  from  p at 

o 

1 2 

a = C)°  to  p + yr  pU  at  a = 90°  and  caused  bv  stagnation  of  the  flow 
o 2 o - 

at  this  point; 

2)  that  there  is  a high  pressure  region  spanning  nearly  160°  of  arc  along- 
the  lower  half  of  the  circular  cross  section  in  contact  with  a flat  sur- 
face ; 


3) 


that  there  is 
of  arc  along 


a relatively  deep  low 
the  top  of  the  forearm 


pressure  region  spanning 
such  that  theie  results 


nearly  200° 
a rather 


# 
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Figure  8.  Surface  pressure  distribution  around  forearm  having  two  degree  angle 
and  resting,  but  not  pressing  against  an  arm  rest. 


substantial  net  force  that  acts  to  lift  the  cone  (or  dislodge  tlie  limb) 
from  tiie  surface  with  wliich  it  is  in  contact; 

4)  tiiat  pressure  minima  occur  near  tlie  top  of  tiie  limb  and  move  lee-ward 
with  decreasing  angle  of  attack  at  a rate  mildly  dependent  on  the  cone 
half-angle.  These  pressure  minima  correspond  roughly  to  points  near 
where  the  flow  will  separate,  and  their  value  lies  between  p and 

5 2 ” 

Po  - PlJy  > depending  on  a. 


Slid  ION  VII 

THE  EFFECTS  OF  FLUID  VISCOSITY 


The  significance  of  fluid  viscosity  in  the  real  world  can  be  dramati- 
cally illustrated  by  a simple  calculation.  Let  us  consider,  for  example, 
the  total  cross-force  per  unit  length  of  cone  predicted  by  the  inviscid 
theory  of  Section  VI,  in  the  special  case  n = 1.  For  this  case,  the  cross- 
force per  unit  length  may  be  written: 
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where 


Po 


surface 


is  defined  by  equations  (52)  and,  more  generally,  (54),  i.e. 
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surface 
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(82) 
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4 sin  0+2  cos  6 tan  3 sin  2a  + C 


Pa=0 


(83) 


Substituting  equations  (82)  and  (83)  into  (81)  and  performing  the  indicated 
integration  leads  to  the  result: 
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R sin 


2ct  tan  6 


(84) 


which  is  a force  acting  in  the  x-y  plane,  located  midway  between  the  normal 
to  the  axis  of  revolution  of  the  cone  and  the  normal  to  the  wind  direction 
(see  Allen  and  Perkins,  1951).  Examining  equation  (84)  more  closely,  one 
sees  that  the  inviscid  theory  predicts  that  f = 0,  i.e.,  there  is  rm  net 
cross-force  acting  on  the  limb  when  8=0°  (that  is,  if  the  truncated  cone 
reduces  to  a simply  right-circular  cylinder)  or  when  a = 90°  (that  is,  if 
the  truncated  cone  has  its  axis  at  right  angles  to  the  incident  stream). 
Clearly,  this  result  is  nonsense,  for  it  says,  in  effect,  that  there  is  no 
such  thing  as  drag.  This  type  of  reasoning,  together  with  some  simple  ex- 
perimental observations,  is  what  led  Prandtl,  at  the  turn  of  the  century. 


37 


to  formulate  his  classic  theory  of  the  boundary  layer.  Essentially,  Prandtl 
argued  that  no  matter  how  small  the  viscosity  of  a fluid  was,  and  no  matter 
how  high  the  Reynolds  number  of  the  flow  was,  the  effect  of  fluid  viscosity 
could  not  be  neglected — that  is,  tiie  flow  could  not  be  considered  to  be 
inviscid  in  the  immediate  vicinity  of  a bounding  surface.  Physically,  this 
is  due  to  the  nature  of  frictional  interactions  which  take  place  between 
fluid  and  surface,  leading  to  significant  losses  of  kinetic  energy.  The 
lost  energy  becomes  critical  if  the  fluid  flows  against  an  adverse  pressure 
gradient,  i.e.,  if  the  pressure  increases  in  the  direction  of  flow,  as  it 
does  on  the  lee  side  of  blunt  bodies  (see  Figures  b,  7,  8 and  9 for  9 > 

90“ ).  In  such  cases,  the  fluid  nearest  to  the  surface  may  have  insufficient 
energy  to  travel  "uphill"  and  it  therefore  either  conies  to  rest  or  reverses 
its  direction  of  motion  under  the  influence  of  the  adverse  pressure  grad- 
ient. When  this  happens,  the  fluid  still  moving  downstream  losses  contact 
with  the  bounding  surface  and  we  say  that  the  flow  has  "separated",  or 
failed  to  follow  faithfully  the  contour  of  the  boundary  to  which  it  was 
originally  attached. 

The  phenomenon  of  flow  separation  is  the  most  important  aspect  of  high- 
speed flow  over  blunt  bodies  of  revolution.  Unfortunately,  it  is  also  the 
most  difficult  to  handle  analytically  because  of  the  nonlinear  elliptic 
nature  of  the  problem.  Thus,  much  of  the  work  in  this  area  during  the  pre- 
large-high-speed-computer  years  was  experimental,  and  the  theories  formu- 
lated were  semi-empirical.  For  instance,  to  correct  equation  (84)  for 
the  observed  effects  of  fluid  viscosity,  an  additional  term  was  simply  added 
linearly  to  the  potential  flow  solution.  This  term  included  an  experimen- 
tally determined  "drag  coefficient",  defined  as: 


s = 


"viscous"  drag  force  per 
unit  length  of  surface 
1 ,;2 
2 PV 


per  unit  of  projected  area 
normal  to  the  flow. 


(85) 


In  equation  (85),  the  "viscous"  drag  actually  includes  two  effects  — 
one  is  due  to  the  rubbing  or  shearing  action  attributable  to  friction  be- 
tween surface  and  fluid;  the  other  is  due  to  separation  of  the  flow  itself, 
which  results  in  what  is  called  a profile,  or  form  drag.  Profile  drag  is 
not  a shearing  action.  It  is  a net  drag  force  which  results  from  a high- 
pressure  area  on  the  windward  side  of  a bluff  body  and  a low-pressure  area 
on  the  lee  side,  the  latter  being  a direct  consequence  of  flow  separation. 

In  other  words,  the  high  energy  losses  associated  with  separation  of  the 
flow  do  not  allow  for  the  complete  pressure  recovery  shown  ideally  in 
Figures  6 through  9 by  the  symmetry  of  the  cur  es  with  respect  to  their 
respective  minimum  value.  Instead,  there  is  a significant  loss  of  pressure 
energy  on  the  back  side  of  a blunt  body,  such  that  the  pressure  distribution 
in  the  separated  region  of  flow  no  longer  balances  that  on  the  front  side  of 
the  body,  and  there  is  a net  pressure  drag  due  to  viscosity.  In  the  problem 
being  considered  here,  profile,  or  form  drag  is  by  far  tlie  dominant  factor 
in  equation  (85).  This  was  discussed  earlier  in  Section  111  (Preliminary 
Assumptions)  and  is  important  in  that  it  allows  to  be  handled  analyt- 
ically as  described  later . 
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Recalling  that  V = U sin  a,  and  noting  that  the  projected  area  normal  I 
to  the  tlow  per  unit  lengtli  ol  a cone  of  radius  K(x)  Is,  locally,  2R(,x)  , we  1 
may  write  tlie  viscous  drag  force  in  equation  (85)  as:  I 


= pU  ^R(x)C  (x)sin^  a 
L)  o 1) 


per  unit  length 
of  surface. 


(86) 


where  ^|j(^)  is  the  "local"  drag  coefficient  at  station  x,  meaning  that  it  is 

a function  of  R(x) , as  well  as  the  local  mean  Reynolds  number,  pVR(x)/,.,  and 
the  Mach  number,  V/V*.  V*  represents  the  local  speed  ot  sound  and  u is  the 

viscosity  of  the  fluid.  Observe  that  tills  force  is  in  tlie  direction  of  V, 

such  that  the  tendency  is  for  the  flow  to  pull  the  blunt  object  alcug  with 
it.  Bearing  in  mind  the  discussion  of  Section  VI  relative  to  an  object  in 
contact  with  a flat  solid  surface  (where  it  was  shown  tliat  there  is  a net 
force  tending  to  lift  the  object  from  the  surface),  we  may  now  conclude 

that  the  forearm  In  contact  witli  a seat  arm  rest  is  subjected  to  forces 

which  tend  to  pull  it  up  (away  from  its  support) and  out  (away  from  tiie 
thorax)  due,  respectively  to  the  generation  of  lift  and  drag  forces  in  tlie 
configuration  described.  This  result  is  in  agreement  with  wind  tunnel  ex- 
periments, which  liave  shown  a distinct  upward  and  outward  motion  of  the 
forearm  when  subjected  to  conditions  similar  to  those  simulated  liere 
(Brinkley,  J.  W. , personal  communication,  and  1974  reference). 


The  Determination  of 


In  general,  tlie  drag  coefficient  is  an  experimentally  determined  quan- 
tity. Allen  and  Perkins  (1951)  suggest  a means  for  estimating  this  quantity 
for  a slender  body  with  varying  cross-sectional  area.  They  first  define  an 

effective  diameter,  D'  = — , where  , the  plan-form  area  is  given  by 


2 


L 

R(x)dx,  and  L is  the  total  length  of  the  body. 


Then,  they  obtain  a 


o 


tabulated  drag  coefficient  for  a circular  cylinder  of  diameter  D'  at  the 


U sin  a(pD') 

cross-flow  Reynolds  Number,  , and  the  cross-flow  Mach  Number, 

U sin  a 

^ . Finally,  they  correct  for  the  effect  of  the  finite  fineness 

ratio,  as  discussed  (and  for  which  tabulated  data  is  presented)  by 

Goldstein  (1938,  pp.  418-421,  and  439)  and  assume  that  the  slender  body 
behaves  with  essentially  the  same  characteristics  as  this  fictitious  circu- 
lar cylinder  of  constant  diameter  I)'. 


I 
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Kflly  (1954)  relined  the  approach  of  Allen  and  I’erkins  to  include 
transient  effects  related  to  flow  development  over  the  blunt  body.  Ho  de- 

fines  a dimensionless  parameter,  — t = — = tt  tan  u,  such 

R R U cos  u R 

o 

that  : 


m=«>  2m+l 

X . 2m+l 

C = S C --  — tan  a 

D m „2m+l 

m=o  R 


where  tlie  various  C are  still  experimentally  evaluated  coefficients, 
m 

In  1959,  Mello  (1959)  pointed  out  that  the  results  of  Allen,  Perkins 
and  Kelly  were  all  strongly  dependent  on  the  body  fineness  ratio  and  the 
angle  of  attack,  due  to  their  basic  inability  to  adequately  predict  the 
magnitude  and  distribution  of  forces  resulting  from  flow  separation. 
Following  reasoning  similar  to  that  employed  in  justifying  assumption  2 of 
section  III,  Mello  therefore  proposed  to  model  theoretically  the  separated 
flow  regime  on  the  leeward  side  of  bodies  of  revolution.  His  model  con- 
sisted of  simple  line  vortices  symmetrically  situated  (along  the  z-axis  of 
Figure  2)  on  the  lee  side  of  a circular  cylinder  in  a uniform  stream,  and  he 
utilized  potential  flow  theory  to  define  mathematically  the  corresponding 
flow  field.  Mello  was  tlien  able  to  show  that  a simple  relationship  existed 
between  the  strength  of  the  vortices  and  the  viscous  normal  force  acting  on 
the  cylinder,  i.e., 

r 

2ttU  c 
o 


where : 

r = the  vortex  strength, 

U = the  free-stream  velocity, 
o 

c = the  radius  of  the  cylinder, 

C,,  = the  viscous  normal  force  coefficient  for  all  the  viscous 
N 

^ normal  force  behind  station  x (see  Figure  2), 

z = z-coordinate  of  right-hand  real  vortex, 

= z-coordinate  of  right  hand  image  vortex, 

the  subscript  x refers  to  conditions  at  the  axial  station  x,  and  "normal" 
means  perpendicular  to  the  cylinder  axis,  x (i.e.,  in  the  direction  of  V). 
Equation  (88)  thus  defines  as  the  drag  coefficient  due  to  profile,  form 

or  pressure  drag,  which,  as  previously  described,  may  be  used  to  approximate 
tlie  cross-flow  drag  coefficient  due  to  viscosity  in  equation  (86).  However, 
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ill  order  to  determine  one  must  know  the  striMigth,  ; , and  loention  (z  - 

z^)^  of  vortices  bcliind  a cylindrical  body.  Again,  much  experimental  work 

has  been  done  in  this  area  and  analytic  studies  have  begun  to  surfate  ,i- 
well.  Notable  among  these  is  tlie  mathematical  formulation  of  Sarpkava 
(1968),  who  lias  attempted  to  describe  the  physical  mechanism  wliich  controls 
the  separation  and  growth  vortices  behind  a bluff  body  by  the  use  of  a 
sophisticated  potential  flow  model.  Sarpkaya  presents  his  results  in  terms 
of  a parameter  k,  where  k lies  between  1.0  and  2.0: 

r = K At  (89) 

The  characteristic  time  interval,  At,  is  obtained  from  a consideration  of 
the  variation  of  pressure  across  the  boundary  layer  at  the  point  of  flow 
separation,  and  V = before.  V'ortex  distribution  behind  the 

cylinder  for  a flow  started  impulsively  from  rest  was  calculated  using  a 
numerical  iterative  scheme,  and  early  results  compared  favorably  with  the 
experimental  data  available  at  the  time.  Such  data,  however,  is  seriously 
lacking  for  the  type  of  geometrical  configurations  depicted  in  Figure  1, 
and  this  would  therefore  suggest  an  important  avenue  of  research  to  be 
pursued  in  future  AMRL  efforts. 

In  any  case,  it  appears  to  be  quite  feasible  to  extend  the  analytic 
work  developed  in  this  report  to  include  the  viscous  drag  force  which 
results  from  flow  separation.  Vortices  of  known  (or  assumed)  strength  may 
be  placed  on  the  lee  side  of  any  of  the  configurations  thus  far  examined, 
and  their  respective  potential  functions  added  linearly  to  equation  (28) 
for  the  uniform  stream.  The  surface  pressure  distribution  may  then  be 
readily  calculated  from  Bernoulli's  equation  once  the  velocity  field  given 
by  the  complex  velocity  potential  is  known.  And  once  the  pressure  distribu- 
tion is  known,  the  forces  tending  to  displace  the  limb  from  its  ejection 
position  may  be  determined  in  a straightforward  manner  — either  purely  by 
analytical  means  or  by  numerical  methods  if  the  mathematics  gets  too  compli- 
cated . 


SECTION  VIII 
CONCLUDING  REMARKS 


It  has  been  shown  in  this  report  that  one  can  approach  the  general 
problem  of  limb  flailing  during  ejection  by  constructing  a simplified 
model  from  which  the  aerodynamic  forces  exerted  on  an  articulated  body  sub- 
jected to  windblast  may  be  reasonably  estimated.  This  model  can  be  t rogres- 
sively  refined,  made  more  sophisticated,  and  streamlined  to  provide  impor- 
tant data  for  both  experimental  and  numerical  studies  of  human  body  dynanilcs. 
In  the  early  form  developed  here,  the  body  joints  themselves  have  not  been 
defined  and  it  would  be  desirable  to  look  into  tliis  in  greater  detail  in 
order  to  assess  the  significance  of  finite  length  segments,  end  effects, 
three-dimensional  contours  and  interactions  wliicti  take  place  between  the 
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limbs  us  a I'csulL  of  their  joint  connections.  Along  the  same  lines,  the 
velocity  U tiuist  be  viewed  as  the  resultant  of  an  outer  flow  that  h.'s  under- 


gone specific  interactions  with  other  body  segments  before  it  actually  im- 
pinges upon  a given  limb,  and  this  is  an  area  wliich  has  not  yet  been 
explored.  We  mentioned  it  briefly  at  the  end  of  Section  II  with  respect  to 
equation  (11)  and  again  after  equation  (55)  in  Section  VI,  when  the  presence 
of  a torso  was  discussed  relative  to  the  angle  of  attack  at  which  the  fluid 
approaches  the  forearm. 


In  refining  the  tlieory,  one  might  also  back-track  to  examine  more 
closely  the  consequences  of  other  assumptions  made  in  Section  111.  For 
example,  the  roughness  of  the  surface  over  which  the  fluid  flows  (which  niay , 
in  part,  be  due  to  the  type  of  ciothing  being  worn  by  the  pilot)  may  lead  to 
certain  types  of  turbulent  behavior  that  is  not  negligible.  Similarly,  at 
.Mach  numbers  exceeding  0.5  compressibility  effects  should  certainly  be 
included  at  some  point  in  the  analysis.  And  finally,  although  the  present 
report  has  ct)ncent rated  specifically  on  the  forearm,  the  theory  could  be 
developed  in  a more  general  form  for  straightforward  extension  to  the  other 
limbs  of  the  body  and  to  other  body  positions  (e.g.,  holdin  a D-ring,  one 
limb  pressing  against  another  - as  the  upper  arm  against  the  torso,  etc.). 
The  foundation  for  doing  this  lies  in  Section  II  (The  Analytic  Body  Model) 
and  in  the  generalized  complex  velocity  potential  of  Section  V. 


For  the  immediate  future,  tiowever,  it  would  seem  that  the  inclusion  of 
flow  separation  effects  as  discussed  in  Section  VII  should  be  the  next  logi- 
cal step  in  the  mathematical  analysis.  Separation  of  the  flow,  together 
with  the  appearance  of  stagnation  points  in  the  fluid  present  themselves  as 
first  order  effects  responsible  for  the  generation  of  limb-dislodging  forces 
that  result  in  flail  injuries.  The  former  gives  rise  to  a pressure  gradient 
that  tends  to  drag  the  limb  along  with  the  air  flow,  and  the  latter  results 
in  a pressure  gradient  that  tends  to  separate  the  limb  from  a surface  with 
which  it  is  in  contact.  Both  of  these  effects  are  functions  of  the  square 
of  tiie  velocity  of  the  incident  stream  (see  equations  (89)  and  (76),  for 
example)  and  so  they  become  rather  significant  at  speeds  near  sonic. 


UST  iiK  SYMWJLS 


Forearm  Vector  from  Klhow  Joint,  A,  to  Wrist  Joint,  IS. 


Plan-Form  Area  = 2 R(x)  dx. 


Total  viscous  dra>4  coefficient. 

I.ocul  Drag  Coefficient  at  station  x. 

Viscous  Normal  Force  Coefficient  for  all  the  viscous  normal  force 
behind  station  x;  Drag  Coefficient  due  to  profile,  form,  or  pres- 
sure drag. 

Bernoulli  Constant. 

experimental  Coefficients  in  the  series  expansion  of  C . 


Surface  Pressure  Coefficient  = 


P ~ Pq 
1 2 • 
2 P^o 


Surface  Pressure  Coefficient  due  to  fluid  of  velocity  U impinging 
upon  conical  surface  lying  parallel  to  the  stream. 

Effective  Diameter  = Aj,/L. 

The  point  z = 0,  y = c. 

The  point  z = 0,  y = -c . 

Chord  line  representing  in  end  view  the  surface  of  intersection 
between  the  forearm  and  seat  rest  against  which  it  is  pressing. 

Total  length  of  arbitrary  body. 

Arbitrary  point  in  space. 

Forearm  cross-sectional  radius  at  location  x . 

o 


Velocity  of  air  stream  approaching  forearm  = U^i  + l'2j  + 
-► 

x-component  of  U,  parallel  to  limb  centerline. 


Magnitude  of  U 
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X ' -cumpiment  ot  L’ , 


y ' -compoiitMiL  ol  U. 


z ' -'omponfnt  of  U. 


y-compoiient  of  U , cross-flow  velocity. 

Local  speed  of  sound. 

Radius  of  circle  tangent  to  y-axis  = a(x)  for  a cone,  (n  = 0). 

Radius  of  free  cylinder  (=R(x)  for  n = 1);  Also,  the  coordinates 
of  the  intersection  of  the  truncated  cone  with  the  y-axis,  which 
constitutes  the  chord  line  GE. 

Cross-force  per  unit  length  of  surface,  not  including  the  effects 
of  flow  separation. 


f,(x) 


Viscous  drag  force  per  unit  length  of  surface. 


►-1  ; Also  used  as  a superscript  (from  1 to  15)  to  designate  a 
particular  body  segment;  Also  used  as  a subscript  to  designate 
that  the  quantity  is  the  image  (e.g.,  a vortex  image)  of  a 
quantity  previously  specified. 

Unit  vector  along  the  direction  of  x'. 

Subscript  for  x-y-z  axes,  i.e.,  j = x,  y,  z. 

Unit  vector  along  tlie  direction  of  y'. 

Subscript  used  to  designate  a particular  body  joint. 

Unit  vector  along  the  direction  of  z'. 

|AB|  = Length  of  forearm  from  elbow  to  wrist. 

Subscript  used  as  an  index  for  the  coefficients  C . 

m 

Used  as  a general  index  for  Integers  0,  1,  2,  ... 

Also  used  as  a parameter  in  the  complex  velocity  potential  of 
Section  V. 

Used  as  a subscript  to  designate  free-stream  conditions 


P 


Fluid  pressure 


Freo-rit roam  pressure 


Surlace  pressure  at  any  axial  station  when  a truncated  cone  con- 
figuration lies  with  its  axis  parallel  to  an  incident  stream  of 
fluid . 

Magnitude  of  fluid  velocity  at  a specified  point  in  the  flow, 

2 2^2 
q = u + V 

2 2 2 

Radial  coordinate  in  cylindrical  coordinate  system,  r = y + z . 

Radius  vector  drawn  from  the  origin  of  an  inertial  coordinate 
system  to  the  k'th  body  joint  (1  = elbow,  A;  2 = wrist,  B) . 

Time  coordinate. 

Magnitude  of  fluid  velocity  in  the  y-direction  at  an  arbitrary 
point  in  the  flow. 

Magnitude  of  fluid  velocity  in  the  z-direction  at  an  arbitrary 
point  in  the  fiow. 

Velocity  vector  at  an  arbitrary  point  in  the  flow  = a vector 
having  components  u and  v along  y and  z,  respectively. 


Complex  velocity  potential  = (j  + iip. 

Complex  Conjugate  of  the  velocity  potential  w;  = cj)  - iilj. 


Coordinate  axis  constituting  the  center-line  of  the  limb  when  it 
is  free,  or  the  line  of  contact  between  the  limb  and  a solid  sur- 
face . 


One  of  the  coordinate  axes  of  an  inertial  coordinate  system. 
A specific  location  along  the  limb. 

Coordinate  axis  chosen  so  that  the  fluid  approaches  the  limb 
parallel  to  the  y-axis,  in  the  x-y  plane. 

A second  coordinate  axis  in  an  inertial  coordinate  system. 

Coordinate  axis  perpendicular  to  the  x-y  plane, 
x-y-z  coordinate  system  is  fixed  to  the  limb. 


The  third  coordinate  axis  in  the  inertial  coordinate  svstem. 
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A 


► 


I 

I 


[ 


An^lt*  ot  ilt.'ick  U relative  ti>  limb:  cos  Jt(t) 


V 


Tan  a = rr 


^orearm  riine  half-anale:  Tan  fi  = 

(Jx 


.] 


Stream  Fu net  ion  of  the  flow. 

Potential  Function  for  the  velocity  of  the  flow:  v = - Vif. 

Unit  vector  defining  the  j'th  coordinate  axis  (j  = x,  y,  z)  of 
the  i'th  body  segment  (i  = 1 to  15)  as  a function  of  time. 


Unit  vector  defining  the  direction  of  the  air  stream  approaching 
the  forearm:  ll/U  . 


i.n  (r^/Cj^),  where  r^  and  r^  are  intersecting  radius  vectors  from 

points  G and  E,  respectively,  defining  the  point  P in  space. 

This  is  one  of  the  coaxal  coordinates. 

A parameter  having  a value  between  1.00  and  2.00  and  acting  as  a 
proportionality  constant  between  the  vortex  strength,  the  square 
of  the  fluid  velocity  and  a characteristic  time  interval. 

Fluid  viscosity. 

Fluid  mass  density. 

Azimuthal  coordinate  in  cylindrical  coordinate  system,  tan  0 = 


The  angle  0 is  measured  counterclockwise  from  the  positive  y- 
axis . 


C ©1  - 02,  where  0^^  and  02  are  the  angles  which  r^^  and  t2  make, 

respectively  with  the  y-axis.  This  is  the  second  of  the  two 
coaxal  coordinates. 


6 

The 

angle  - 

cot  0 . 

X 

The 

complex 

variable  y + 

iz  = ic  cot  Y 

ie 

= re 

X 

The 

compl ex 

conjugate  of 

x;  = y - iz  = 

-ic  cot  j = re 

4 

Tlie 

complex 

variable  C + 

in. 
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n|  >> 


I 

I 


i 


V 


oo 


r 

z 

A 


Tilt?  complex  conjugate  of  i;;  = ^ - iri. 

The  Gradient  operator. 

Designates  infinity. 

Tile  strength  of  a vortex. 

Designates  a summation. 

Designates  a difference  between  tlie  twt>  quantities  following 


r 


IT 


3.14159. 


COMPUTER  PROGR.\M  FOR  CALCUEATINC  AND  PEOTTINC  STREAM  FUNCTIONS 
DEFINED  BY  EQUATION  (ib)  IN  THE  TEXT  MATERIAL 

APPENDIX 


DLB,T350,CM100000,STCSB.  L7305J3  BRUNGART 
ATTACH , D ISSPI.A , 1D=X654  32  1 , SN=AF1T  . 
L1BR.\RY,DISSPLA. 

FTN. 

EGO. 

ONLINE. 


PROGRAM  FLOW ( INPUT .OUTPUT ,PLF I LE=0) 
COMMON  YLAST.XLAST 
CALL  CO.^ffRS 
CALL  BGNPL(l) 

CALL  CROSS 

CALL  TITLEdH  ,-l,"Y",I 
CALL  GRAPH(-5.  ,I  . ,0.  ,1.) 

CALL  CIRCLE 
PI=355/113. 

A=1 

U=1 

API=A*PI 

APU=API*U 

L=2 

DO  6 K=l,2 

DO  6 .M=l,701,25 

CALL  BLNKl(8.9,10.0,0.  ,7.) 

CALL  BLNK2(0. ,1. ,0.  ,7.) 

IF(L.EQ.2)  CALL  BLNK3 (8 . 5 , 9 . 0 , 0 . , 7 . ) 
F=-l 

IF(K.EQ.2)  F=1 
PSIBASE=M*.005 
KK=0 

PSIU=PSIBASE+.OI 
PSIL=PSTBASE-.OI 
DO  5 1=1,200 
Z=. 025*1 
DO  4 J=l,200 
Y=.025*(J-1)*F 
1 RSQ=Y*Y+Z*Z 

F1Y=APT*Y/RSQ 

F1Z=API*Z/RSQ 

F2=TAN(F1Z) 

F3=C0SH(F1Y) 

F4=TANH(F1Y) 

PSI=APU*F2/((F3*F3)*(F4*F4+F2*F2)) 

IF (PST. LT. PST  BASE)  GO  TO  4 


Z", 20, 10.  ,7.) 


000160 

000170 

000180 

000190 

000200 

000210 

000215 


000225 

000226 

000240 

000250 

000260 

000270 

000280 

000290 

000310 

000320 

000330 

000340 

000350 

000360 

000370 

000380 
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IF(((Z-I)*(Z-1)+Y*Y) .LT.l)  GO  TO  4 000390 
IF(PSI  . LT.PS1U.AND.PS!  .GT.PSTL)  CALL  DR.AW(Z  ,Y  , KK)  000400 
IF(PSl.LT.PSIU.ANU.PSl.GT.PSIL)  GO  TO  5 000410 
CONTINUE  000420 
CONTINUE  000430 


CALL  RESKTC'BLNKS'') 

IF(L.EQ.2)  GO  TO  10 
CALL  MIXALFC "GREEK") 

IF(K.EQ.2)  CALL  MESSAG("(Y)  =" , 5 , 9 .0 ,ZLAST) 
1F(K.EQ.2)  CALL  REALNO (PSTBASE , 3 , "ABUT" , "ABUT") 
L=2 

GO  TO  6 
L=1 

CONTINUE 
CALL  ENDPL(O) 

CALL  DONEPL 

STOP 

END 

SUBROUTINE  DRAW(Z,Y,KK) 

COMMON  YLAST , ZLAST 
Y=Y+5  . 

IF(KK.EQ.O)  YLAST=Y 
IF(KK.EQ.O)  ZLAST=Z 
IF(KK.EQ.O)  CALL  STRTPT(Y,Z) 

CALL  CONNPT(Y,Z) 

KK=1 

RETURN 

END 

SUBROUTINE  CIRCLE 
DIMENSION  ZZ(lOl) ,YY(101) 

DO  I 1=1,50 
ZZ(I)=.04*I 
ZZO0PI)=2.-ZZ(I) 


YY(I)=SQRT(1.-(ZZ(I)-1)*(ZZ(1)-D) 

000630 

YY(50+I)=-YY(I) 

000640 

CONTINUE 

000650 

DO  2 1=1,100 

000660 

CALL  C0NNPT(YY(I)+5.  ,ZZ(D) 
RETURN 

000680 

END 

000690 

000440 


000470 

000480 

000490 


000550 

000560 

000570 

000580 

000590 

000600 

000610 

000620 
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